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Abstract 


rmdmg the optniial shapes of the mechanical components is one of the most visited field of the op- 
timization and also this is a difficult task The mam difficulty lies m the method of lepiesentation 
oi shape in an optimization algoiithm Most optimal shape design pioceduies involve pie-fixmg the 
shape's paiticulaily a mathematical form (say through cubic splines or some polynomials) The task of 
\]\c optimization method is the optimization of various parameters governing the mathematical foim oi 
i(l('iitifymg the optimal value of the co-efficients defining the mathematical functions Although in some 
( ases when the pioblem is simple this method woiks but it cannot be geneialised to find the shape of 
any aibitiaiy shape optimization problem The problems, where the information about the shape is 
not known befoie hand, finding the solution becomes very difficult 

To oveicome this difficulty, the representation scheme where the shapes are given by eithei the 
pic'sc'iice 01 the absence of the small material pieces, is used This way, the mateiial shape can be 
K'pic'sented by a stung of binaiy numbers, where each binary variable decodes to the presence oi 
abs(uice of the mateiial at a specified place A smoothing algorithm is used to remove the difficulties 
like shaip coineis and point contact etc The shapes are analysed by using the finite element analysis 
to find out the sti esses and stiains Constant strain triangles are used as the basic elements for finite 
('lenient analysis A hybrid appioach based on the Genetic Algorithms (GA) is used to find out the 
opthiial shapes Foi the single objective problems the objective is to minimize the weight The optimal 
solution obtained by the GA runs, undergoes a hill climbing local search, which ensuies the convergence 
of the solution to the globally optimal solution 

When the problem is solved for more than one objectives, it becomes a multi-objective optimization 
pioblem Here the objective functions are the minimization of the weight and the minimization of the 
(kffiectioii, which aie conflicting in nature Then, a set of Pareto-optimal solutions is obtained as the 
solution to this pioblem, as one single solution cannot be the optimal solution To solve these types of 
pioblems, a number of multi-objective genetic algorithms (MOGA) are given m literature A specific 
multi-objective genetic algorithm, an elitist non-dommated sorting genetic algorithm (NSGA-II) is used 
Each solution obtained by the NSGA-II search undergoes the local search to give an optimal solution 
A weighted sum stiategy is used to convert the multiple objectives into a single objective The weight 
vec tor is calculated by using the fitness values of the solutions in an adaptive way Two strategies are 
adopted for the weight calculation- one the fixed weight strategy, where weights are calculated on the 
basis of initial fitness values, and the second is the continuously updated strategy, where the weights 
aie found by the continuous updating of the fitness values The number of the solutions obtained after 


tho liybiid appioach may be too large for the designer to take into consideration A clustering approach 
IS used to 1 educe the number of solution to the desired number 

A iiiuiibei of test problems like the design of cantilever plate, simple supported plate lioister, etc 
aic' solvcvl Both the single objective as well as multi-objcctivc cases are handled Each problcin is 
solvc'd twice hist -when the weight of the plate itself is not considered for the calculation of the load 
and sc'cond is the case when the load of the designed plate is also accounted for while calculating tlu^ 
load The results obtained by the application of hybrid appioach to all these problems give us ven 
Intercast mg results In a number of cases the results obtained are very obvious and some times some not 
(vis\ to visualize solutions are found to be optimal This makers the approach very efficient m handling 
the' shape optimization problems Most interesting application of this hybrid approach is the design of 
the' bicycle frame where a number of different bicycle frame shapes are obtained m just one simulation 
1 nil 
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Chapter 1 


Introduction 

1 1 Introduction 

Dc'si”!) oi the shapes of the mechanical components is not a new activity In fact, the shape optimi/.at loii 

IS oiK' of file most visited areas of the optimization The most populai method of shape optimization 

au' the use of classical methods Since many classical methods require the function and the deiivativc' 

mtoi Illation Therefore, the shape of the component is decided to be approximated by a mathematical 

\ 

luiK tioii mostly a polynomial m nature, as the mathematical functions can be easily handled by the 
( lassie al iiK'thods The co-efficients of the polynomial and other parameters governing the mathematical 
hnulion considered are obtained using the different optimization methods, to find the optimal shapers 
[2 11 12 16] This methods suffer from the disadvantage of requiring the knowledge of the shape before 
hand lu any arbitrary problem where the loading and support conditions are such that it is difficult 
to giic'ss the shape, the method finds it difficult to get the shape Besides this, the problems where 
the' properties of the material under consideration are not constant over the search domain, the initial 
guess of the shape is not at all easy The shapes represented by the mathematical functions also find it 
difficult to find the holes m the shapes, which may be helpful m some of the cases as the material can 
be sometimes reduces from the areas of low stress Despite these drawbacks, the classical methods aie 
most popular due to the absence of any other approach which does not suffer from these drawbacks 

Here m this work a representation scheme other than the representation of the shapes by the math- 
('uiatical functions is adopted Here the shapes are derived from the rectangular plate, divided into 
a iiumbei of vsmall elements The shapes are represented by the presence or absence of the material 
('lemeiits [4, 5, 13] Thus each shape is presented by a string of binary numbers Each 1 in the binary 
string signifies the presence of the material and 0 corresponds to the void in that position The basic 
shape IS extracted from this string, through a smoothing techmque, ehmmating the string representing 
the' unconnected regions FEM analysis is used to find out the stresses and strains for these shapes 
A hybrid appioach ba£«d on the genetic algorithms (GAs) is used to solve the real world problems of 
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the shape optimization Binary GA is used to find the optimal shape of the component and then a hill 
( limbing local search is used to solve the problems 

riu' piobkmis can be of optimization of single objectne oi of optimization of multiple objectives with 
a I('\\ (onstiamts like geometiy constiamt and stiess and stiain constraints Foi the multi-objectue 
optiim/ation piobleiiis, the solution is not a single point instead it is the set of solutions whidi die 
solutions known as the Paieto-optimal solutions To find the Paieto-optimal set a iiumbei of 
limit i-()b]('( tive GAs (MOGA) are suggested [6, 15 21 24 29 32, 38] A specific MOGA- -ehtist non 
(louunatecl genetic algoiithm (NSGA-II) is used to find the Pareto-optimal solution set Each solution 
(){ th(^ solution set obtained fiom the MOGA simulations are then improved by the local seaich The 
nsc' of th(' local seaich ensures a better conveigence to the global Pareto-optimal front and also help 
in K'diuing the size of the non dominated solutions set to a reasonable number Finally a clusteimg 
algoi it hill IS used to find the desired number of different solutions The efficacy of the pioposed method 
IS shown tliiough a number of test problems like — the design of the cantilever plate design of simplv 
suppoit('d plate, design of a hoister, etc The results obtained are very interesting as the approach gives 
SOUK' veiA intuitive results and also it is able to find results which are not so easy to visualize The 
])iop()s('d method gives the solutions that show the trend of the placing of the material m the optimal 
luamun Since the solutions obtained are evolved through the GA which mimic the nature, the solutions 
oblaiiK'd aie very near to the global optimal solution And the same is shown with the help of a numbei 
oi tc'st pioblems 

1.2 Previous work in the field 

Til th(' activities related to design of the mechanical components, shape design is the first step [18] 
SiiK e the development of the numerical techniques during the world wax-II, shape design of mechanic al 
(ompoiieiits has also become very popular [2, 11, 12, 31] Most of the studies are classified into two 
luajoi categories - 

1 Shapes are represented as a mathematical function 

2 Shapes are formed by deleting material from an initially chosen shape 

In the fust approach, the design of the shapes is done by pre-decidmg a mathematical function and then 
hiiding optimal parameter of shapes [16, 33, 35] This approach requires the knowledge of the shapes as 
well as the proper range of the variables is to be decided Continuity and compatibility conditions are to 
!)(' satisfied on the boundary points The second approaxdi is more elegant and used more sparingly [14] 
Tlic rectangular plate is taken as the initial shape and a number of elementary shapes such as circle, 
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square, or rectangle etc are deleted from different places The optimization problem is to find the 
optimum place and the size of these elementary shapes This approach is very simple m principle but 
suflc'is fioiu a iiumbci of problems Often the shapes obtained by the subtiaction can not be lealized 
1 he disc oiitiiiuities in the boundaiy shapes are present Since a number of element aiy shapes form the 
hoiuidaiy of the lesultmg components, decision variables are highly inter -dependent 

A number of researchers [4, 5, 13, 19, 22, 23, 28] have shown that the shapes can be represented 
by the binary material/ void representation The use of genetic algorithm to this type of shape opti- 
nii/ation problems gives very interesting results A number of studies are done for the single objective^ 
optiiiii/ation problems [4 5, 13 19, 22] The researchers have used the Finite Element Analysis to 
analyze the solution and find the stress and strains The representation taken by the icseaiclieis suffcis 
tlu' drawback of not considering the cases of point contacts, which may lead to very high stress Then 
t he results of the study may not remain the optimal The use of the genetic algorithm approach to solve 
some real life problems is also given m [28] where the design of a car bumper is carried out using the 
g(m(dic algorithms The use of single objective GA to solve the bicycle frame design problem is carried 
out m the [23] 

But so far nobody has carried out the design of shape optimization problems for the cases when more 
than one objectives are considered This study presents the design solution with the new representation 
which does not have the problem of sharp corners and the point contact for single and multiple objective 
l)i()blcms And the solutions are solved through a hybrid approach combining the advantage of the GA 
and hill climbing local search 

1.3 Organization 

This study is presented in all six chapters The brief description of all the chapters is helpful m studying 
the work in an organized manner This is given here - 

Chapter 1 is the introduction to the work In this chapter, the difficulties of the classical approaches 
in solving the shape optimization and the brief outline of the proposed approach is discussed The 
work so far done by others is also discussed briefly 

Chapter 2 discusses the method how GA can be apphed to the problem of shape design This chapter 
also discusses the details of the FEM analysis, the 2-D crossover operator, and the method of local 
search for the single objective problems 

Chapter 3 presents the description of the different test problems taken into consideration Then the 
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lesults foi the single objective optimization are presented In this chapter the results of comparison 
of the local search method and the suggested hybrid method is also presented 

Chapter 4 beais the brief intioduction to the terminology used m the multi-objective optiniizatiou 
piobleins and then biiefly NSGA-II algorithm is also piesented The local scaich method and 
(liffeient weiglit calculation strategies are discussed 

Chapter 5 piescnts the results for the test problems when the numbei of objectives is moic than on(‘ 
The lesultant Paieto-optimal set and the shapes obtained are piesented In the end of the chaptei 
tli(' 1 ('suits for the design of bicycle frame is presented, which are found to be the most exciting 

Chapter 6 sunimaiizes the hybiid method, mam results and the conclusions drawn fioni the results 
Tlu' (hajitei also gives some of the future scope of the work, which will be useful foi mteiested 
leadeis 

1.4 Closure 

Th(" chapter gives an overview of the study This discusses the drawbacks of the classical methods used 
toi solving the shape design pioblems The scheme of representing the shapes by the binary strings is 
given 111 brief and then outlines of the proposed hybrid method to solve the shape design problems are 
sketcJied foi both the cases of single objective problems and the multi- objective problems The previous 
work done is mostly done for optimizing the shapes formed by the group of rectangular elements, without 
any consideration to the manufacturing problems Most of the work is done for the single objective 
optimization Some interesting practical applications are also reported in history The organization of 
the report helps in guiding the reader to the section of his/her interest 
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Chapter 2 

Genetic Algorithms for single-objective 
optimization 


1 li(^ GouetK Algol ithms (GAs) is a leliable approach of solving the optimization pioblem It niiiiiKs 
t li(' iiat nial ('volutiou process to find the optima [20] It has the potential of finding the optimal solution 
\\!ii(h nici's be suiiouiided by a lot of local optima or it may be highly multi-modal It can be used 
(oi oiluu piiiposes as well, like design of components, scheduling routing etc Here the discussion is 
( ('lit, cued aiouiid the question — How GAs can be applied to the shape optimization pioblems or the 
piolilems of stiuctuies which seem to have no relation with the domain of genetic algorithms Answer 
to this (luostion tells an approach how the problems from different walks of life, can be modelled to 
<t})l)ly G A. This chapter mainly deals with the same aspect 

2.1 standard genetic algorithm 

GeiiotK Algol ithm(GA) is the optimization technique, which mimics the nature m Older to find the 
global optima of a given problem This is inspired firom the concept of "Survival of Fittest” It means 
t hat only those individuals, which are better than others m one or other aspect survive Unlike clcis- 
sH al oplunization methods GA deal with a population of individuals In Binary GA individuals aie 
u'pieseiited by the bmaiy strings called chromosomes These chromosomes are the string of O’s and I’s 
( ailed allele In leal coded GA, chromosomes comprise of string of real numbers Salient features of a 
staiidaid simple GA are as following 

1 Selection 

2 Ciossover 

3 Mutation 
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2 11 Selection 


Vox a gnen population of individuals only those individuals are selected which are bettei in fitness than 
( ){ h(M s Tins way, GA follows the Darwinian law of survival of the fittest and by this selec tion oiii seau h 
IS guided towaids the optima Moie importantly, this method does not reject the solution which is not 
[\\c Ix'st This solution may have the information, useful m finding the global best solution Most used 
iiK'thods foi selection aie 

1 Touinamoiit selection 

2 Stochastic imivei sal sampling 
i Roulette wheel selection 

2 12 Crossover 

III c lossovci two individuals mate and produce children This is a very impoitant opeiator which cieate 
I !i(' new solutions, by using the piopeities of the parents This creates the solutions m the legion othei 
than the existing legion This explores the search space to find the best region It is supposed that the 
individuals which aie suivivmg, have some good string, which let them alive So good parents combine 
to piociuc.e bettei children Two types of crossovers are mostly used 

1 Binaiy crossover 

Deals with chiomosomes with binary numbers 

2 Real ciossovei 

Deals with chiomosomes with real numbers 
2-1.3 Mutation 

This mimics the sudden changes occurring in the nature, that creates the individuals, different from the 
individuals of its own category This operator is particularly helpful when the search has stuck at some 
local optima or search is going slowly This can also create solutions which are not good They will be 
taken care by the selection and other operators 

The algorithm of simple GA is shown in Figure 2 1 


2.2 Representation scheme 

Foi all the different problems a rectangular plate is taken as the basic shape The plate material has 
the physical significance whereas the representation in form of the O’s and I’s do not have any physical 
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Figure 2 1 Simple genetic algorithm 


sigiiihc ance but these strings of O’s and I’s can be processed by binary GAs Thus, heie is a discrepancy 
111 the two representatrons Thrs is avoided by finding the compatibility between these representations 
Foi this puipose the rectangular plate is divided in small regular pieces which are called elements These 
('k'liK'iits ate oui basic building blocks Each element can either be present or absent If the element ma- 
ti'iial is picseiit a value ’1’ is assigned to it and if element is void it is specified as ’0’ This way, material is 
jn ('sc'iited in the form of a sequence of O’s and I’s These shapes in form of the string are processed by the 
GA’s Since the body taken is two dimensional so this string is converted in the two dimensional airay 
by filling the array positions by the string members sequentially from the left as shown m the Figure 2 2 

This is more clear when a material shape is represented m the form of an array and this array is 
( oiiveited into a binary one-dimensional string and vice-versa is also true 

Example: The representation scheme is more clear by the example Here is a randomly created string 
of length 100 


1111000111 0100010111 1001101000 0011110010 1110000011 
1100100000 1010101000 0000100110 0011101110 0001010000 







Two dimensional representation 

rii(' ( onospouding two-dimensional array of size 10 x 10 is given in the Figure 2 3 
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Figure 2 3 Array representation Figure 2 4 Material representation of the 

string 


Material Representation 

The material representation for this 2-D array is given in the Figure 2 4 If the material shape is given, 
one can find what will be the corresponding string and if a string is given the material presentation for 
llu' corresponding string can be easily obtained 






2 2 1 Connectivity and smoothing 

Convc'iting the binary string into the material shape, gives a skeleton of the shape but this shape- 
IS iiol uu'aiiiiiglul unless it is connected Unconnected legion will not be able to take any load and 
hem e without c oiitiibuting any thing significant, it will increase the weight Besides this it will not tie 
,i( ( ouiiti'd while pioducing the object physically In ordei to find the connectivity the biggest coiiiiec ted 
K'gioii (( lustei) IS found To find the largest connected legion connectivity between elements is defined 
using lollowiiig definition 


Definition 1 To define connectivity two terms neighbour and group are defined as following 

1 Neighbour. The material elements which are sharing at least one corner are called neighbouis 

2 Group The collection of all the material elements, connected to each other by means of neigh- 
houi', constitute a group The group is often called as cluster 

The most important property of the group is that if one starts from any of the element 
of the group one can reach any other element of the same group 

2 2 2 Steps for finding connectivity 

C'oiiiu'i ti\ ity foi any element is obtained m following steps 

• One single element for which, material is present and which does not have any neighbour, is teimed 
as a sepal ate group 

• The material element, which have at least one neighbouring material element, belongs to the 
common group same as the neighbour and the group number is smaller of the two group numbers 

• Classify all the elements as the members of the different groups or form different groups, such 
that each element is exactly the member of one group only 

• The gioup which has maximum numbers of elements among all these groups, is found and that is 
taken as our basic shape for that string i e the string gives the shape as identified by the biggest 
group 

The biggest cluster is taJken as the unsmoothed shape and the string is now repaired by substituting 
all the elements which do not contribute to the biggest cluster by ’0’ Thus lamarcfcian approach followed 
here [25] The shape after clustering for the string taken for the example is shown m the Figure 2 6 
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2 2 3 Smoothing the component 

\ft(H oxtiactmg the basic shape from the string a region, connected by the edges oi have coiiiei roii~ 
ta( ts 1 (' the two elements aie connected by just sharing one corner, may be obtained These tvpc^s ol 
( oiifigmation aie feasible foi theoretical works, but for designing of the object undci action of actual 
load IS not possible as the design will be highly unsafe at the point, where the legion is just shaimg onc' 
node' ouh because there is not enough material to support the load and hence the stress may exceed 
t ii(' allow c'd limit So with this representation scheme the shapes with ciossed-ribs sort of allallgement^ 
w liK h aie oftem found to be useful in a number of cases cannot be found 

111 01 del to overcome this problem, smoothing operations are applied Smoothing is done by follow mg 
tlu^ givcm conditions One typical cause of each smoothing method is shown in Figure 2 5 

Condition 1 When there is a hole made by the single element, it is filled up by replacing the void b\ 
a material element 

Condition 2 When a void is having two material elements adjacent to it, it forms a sh^rp corner and 
this will lead to high stresses To avoid it a right angled triangle is put there such that the shaip 
( 01 iiei gets eliminated This filling of the triangle is m line with the concept of placing the fillet 
iKVu the sharp corner 

The use of triangles can be justified by the fact that the fillet can be modelled as a polynomial of 
high degree where as the triangle can be taken as a fillet with a polynomial of degree one 

Condition 3* When the void element is cordoned by three material elements, the triangles are placed 
such that 3/4t/i of the void is filled Which l/Ath part is to remain void is decided by the side of 
th(' vacant element If right hand element is absent, the l/4th portion of the right side can be left 
as void 

Condition 4 When the arrangements like shown in first part of the Figure 2 5 (sharing of one node 
only) are encountered, two right angled triangles are placed such that, a crossed rib like shape is 
obtained 

Aftei smoothing the final shape represented by a string is obtained This shape becomes a candidate 
foi the consideration of finite element analysis, provided it satisfies the geometry constraints, if any 
Tlie final shape represented by the string taken m example is shown in Figure 2 7 
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Figure 2 6 Biggest cluster 


Figure 2 7 Smoothed shape 


2.3 Single objective problems 

Single objective problems are the most basic optimization problems Most of the classical methods 
aie applicable to single objective problems Many classical methods for multiple objectives convert the 
jiioblem into a single objective problem through a number of means and then solve it. So it is essential 






U) find whether this strategy of finding the shapes using Genetic Algorithms based hybrid method works 
w('ll oi not The problem is formulated with an objective of minimization of the weight Constraints 
whidi all the feasible solutions must satisfy, are imposed First constraint is that the stiuctiue must 
not \ lolate the stiess constiaint and second, the maximum deflection must be less than a allowed le\el 

2 3 1 Objective function 

ho] Single objective problems the objective function is the weight of the object, in most of the cases 

Weight = Area * Thickness * Density (2 1) 


1 oi all the practical designs the thickness and the density of the material remains constant foi all the 
individuals and over the generations So the area of the individual shape can be used as a measure of 
( h(' weight Thus the problem requires the weight or the area of the shape to be minimized 

2 3 2 Constraints 

I h(' sha])C' used for carrying the load, must not fail while doing so To avoid failure, the limits on certain 
paiainc'teis are set and these constitute constraints Following constraints, are taken into consideration 

• Geometry constraint Any structure is feasible, only if it satisfies the geometry constraints 
The shape under consideration must have the nodes for applying the load and also it must have 
i.h(' essential number of support nodes If any of these conditions is not satisfied the shape is 
ie]e( ted and a high erior value is assigned 

• Stress constraint Any structure is safe, so far, the stresses arising due to the loading does not 
exceed the maximum stress, the material can sustain The maximum stress arising at any point 
must be less than the allowed stress 

[^max] ^ [^allowed] (2 2) 

In most of the cases allowed stress is same as yield stress Then 

max] ^ [o’yteid] (2 3) 

• Displacement constraint; This constraint is important as the design of shapes is considered 
It has been observed that some structures are safe from the stress point of view but they fail due 
to excessive deflection So it is necessary to keep the maximum deflection m the structure undei 
a specified limit to avoid it from failing due to excessive deflection This is given as 

[^maa:] ^ [^oli owed] (2 4) 
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2 3 3 Material selection 


The piopeities of the material over the elements can be controlled, making it easy to handle different 
n}( ff//s 01 the uai Hit ion in the piopeities This is an advantage of using the small element ( oncept Thus 
ilu' design ot the optimal shapes when the properties of the material aie varying ovei the domain can 
!)(' (aiiK'd out easily, which is otherwise a very difficult job to design By this method diffeient loading 
( oiiditioiis like distiibuted load, axial load or vertical load etc can also be handled easily 

2 4 Fitness and error evaluation 

C.!('iu’tu algoiithm woiks with d population of individuals Each individual is assigned a htness accoiduig 
to its fiuiction value(&) GA tries to maximize the fitness Since this is a constrained optimization 
piobleiii so foi each individual according to its constraint violation the erroi value is assigned Tlu' 
('1101 IS derided by finding the geometry constraint and the stress and the displacement constraints 

2 4 1 Geometric constraints 

The shape obtained after smoothing is not always feasible If the individual does not .have the node 
whole the load is applied, the material will be unstressed, and the finite element analysis to find th(' 
stress 01 displacement will have no meaning If the individual does not have the support nodes, the 
inateiial under t,he action of loads will move Then the problem will be of the dynamics which is not 
^olved here In these conditions the geometry constraint is infeasible and the individual is assigned a 
laigc' eiioi The individual is feasible only if all the nodes where body is supported and the nodes where' 
load IS applied arc present 

2.4 2 Stress and displacement evaluation 

The stresses and displacements are required for the evaluation of the constraints Finite element analysis 
IS used for evaluation of stresses and displacements The details of finite element analysis, are given m 
separate section The maximum stress and maximum strain values at any point appearing m the shape 
under consideration, are checked to find the violation of the constraints 

2.4 3 Fitness evaluation 

The objective is to minimize the weight More is the weight, less is the fitness Since the density of 
the material and the thickness of the plate are taken constant over the plate, the area of the shape can 
be directly taken as a measure of weight The fitness is calculated by finding the plate area of the shape 

Fitness = Maximum plate area — Area of individual shape 
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GA will maximize the fitness 


2 4.4 Constraint handling 

1 li( sK'ps h)i (oustiaints handling aie given as following 

1 S( almg tli(' stic'ss and displacement constiamts 

2 Find tli(' scpiaie sum of the constraint violation (if any) If none of the constraints is violated 
( oiiespondnig violation is set to zero Mathematically the lepresentation is as following 


CVi 


|o~mat| ^alloued 
^allowed 


if 


CTn 


^ \^allowed\ 


NULL else 


(2 5) 


CV2 


’^‘^s2iowed‘’“^ ^ ^allowed 

NULL else 


(2 6 ) 


\\ Ik'U' 


CVi = Constraint Violation 


wlu'K' 


k = number of constraints 

3 Tom nament selection is done to handle the constraints as suggested by Deb, 1999, [9] Tournament 
size taken is two Winner of the tournament between two solutions is decided by following the 
( oust! amt domination conditions 

Condition 1* If any of the individual is feasible and other is not, select the feasible individual 

Condition 2. If both individual are infeasible select the individual with least error or constraint 
violation 

Condition 3: If both individuals are feasible or both are infeasible with same constraint violation, 
select the individual with better fitness 


This IS shown elsewhere [9] that this strategy can handle the constraints efficiently 


2 5 Types of operators used 

Foi (liffcTcut opeiations of genetics following operators are used 

• Tournament Selection is used foi selection Tournament size of two is set 

• Binary 2-D Crossover is devised and used The details of this ciossovei aie given in sepaiate 
s(H tion 

• Siinplo landoin mutation opeiatoi is used heie 

2 6 Hybrid method for shape optimization 

( t \ s (l(v\l with the population of individuals It has the potential to reach the global optimal solution 
hut foi this the population size must be sufficient [17] In these problems each feasible individual has 
lo iindmgo the finite element analysis, which is computationally very costly So the population size is 
takcm small to leduce the computation time The reduced population size may cause conveigence to the 
!()( al optimal solution m a limited number of generations or it may prolong the seaich To leduce the 
nuK' to get the best solution and to ensure the convergence of the solution to global optimal solution 
ilu' knowledge of schema processing by GA’s is exploited Schema processing means that good strings 
< oiiK' closei to find the optimal solution GA’s do schema processing to bring the shoit good schema’s 
logc't hei 

The classical methods, work very well when the search space has the ummodal optima oi if the 
soluf ion starts with a veiy good initial solution This will not work very good and has a laige probability 
of gc'ttmg stuck in the local optima, if the search space has too many optima To ensure the conveigence 
I o the global optima both of these concepts are applied here Initially, the GA’s start with the randomly 
genoiated population and process it to get the good arrangement of the strings This takes the search 
lo some optimal solution, which may not be the global optima Now, the hill climbing local search 
method similar to the steepest descent method is apphed, with the solution obtained fiom the GA as 
1 hc' initial solution This way it starts with a good initial solution and taJees the search to the globally 
optimal solution quickly The combined strategy is given following 

1 Genetic algorithm is used to find the best solution. 

2 Further refinement using one bit hill climbing strategy This hill climbing is a local search method 

The combined strategy has yielded very good results and proposed as a generic method of finding the 
global optimum method for the shape optimization problems 
Hill Climbing Strategy 
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step 1 Start from the left end of the binary string representing the domain and replace the mutai 
the ciiuent bit i e replace ’0’ by ’1’ and ’1’ by a ’0’ 

Step 2 Find out the biggest cluster and find connectivity 

Step 3 II tli(' solution satisfies the geometry constraint, calculate stress and displacenieut b\ finite' 
(dement analysis Find the weight of the solution 

Step 4 If the weight is improved and all the constraints are satisfied accept the change else lestoic' 
tli(' (liaiige and goto Step 1 unless the end of the string is not leaclmd 

2 7 Finite element modelling (FEM) 

The finite element modelling is the most important and computationally costly part of this whole 
(lesigi) pioccss The finite element analysis is carried out on the shape obtained by the decoding of the 
stung The finite element analysis gives the displacements as the primary variables These lesults au' 
l)ost piocessed to find the stresses developed m the body and the maximum stress developed in the 
l)ody The maximum displacement appearing m the body is also extracted from the data Here 

• Maximum stress is obtained by using the Von-Mises Criterion 

• Displacement is obtained by finding the vector sum of the displacements in the x and y diiections 
at ea( h node Maximum displacement is the maximum of these displacements 

• The mateiial is assumed to be ductile 

• Constant stiain triangle is used as the basic element for the finite element modelling 

2 7 1 Steps in finte element model formulation 

The finite element analysis is based on the Principle of Virtual Work This states 

If a general btructure, which is in equilibrium with its applied forces, is subjected to a set of small, 

compatible mtual displacements, the virtual work done by the external forces is equal to the virtual 

bfrain energy of internal stresses 

Mathematically it can be given as 

5Ue = SWe (2 7) 

wheie 

5Ue virtual strain energy 

SWe virtual work done by external forces 

For the finite element analysis the steps given below are followed 
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EleiTient formulation 


riu' two dniieiisional legion shown in Figure 2 8 is divided into straight-sided triangles Foi this eadi 
mtc'iioi luateiial element is divided iii two tiiangles and those lectangulai elenieiits, wliidi aic' on the 
houiidaiv (iiu hiding the boundary of the hole) aie divided into foui triangular elements as shown m 
l-hguK' 2 9 The triangles are also divided into srnallei triangular elements Nodes are placed on all tlu' 
v(>i lev points 

I hi s/.e of tiiangips need not be same The boundai> can be more accurately modelled by putting 



Figure 2 8 Rectangular domain Figure 2 9 Meshed domain 

mou' elements or a refined mesh there The elements of the different sizes are present in the same mesh 
Shape function and degrees of freedom 

Constant strain triangle elements are used for finite element modelling For this the shape functions aio 
lineal over the element The iso-parametnc representation and shape functions m natural co-ordinates 
aio used This gives the freedom to use different sized elements while finite element analysis 

• One dof is assigned to displacement m each direction 

• This IS a two dimensional problem so each node has two degrees of freedom One dof is assigned 
to the displacement m each direction 

Global node numbering and connectivity matrix 

Attei creating the elements the global node numbering is done of all the nodes Global node numbering 
IS done by tracing all the nodes from the top left corner to right bottom comer as shown in Figure 2 10 
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Figuie 2 10 Global node numbering 


node nurnbeiing is done, connectivity is sought i e finding which element is attached to which ele- 
iiumt Foi this piupose moving from the top left material element, connectivity is found and accordingly 
the Connectivity Mafnx is formulated 

Assigning the geometric and material properties 

Tlie piopeities of the domain aie assigned to each element For geometiic properties, the co-ordinates 
lufoimation is stoied and for mateiial properties Young’s modulus and density can be assigned to each 
individual element This helps in designing, when the search domain has more than one elements or 
the mateiials have variation m properties over the domain 

Assembly 

The difierential equations are formulated and solved to get the matrix form By this local element 
stiffness matrix and similarly local element force vector is also obtained Assembly of the local stiffness 
matiices and force vector is done using the connectivity matrix This gives the global stiffness matrix 
and global force vector 


Apply boundary conditions and solve the simultaneous equations 


I 

’ii 

I 





Boundaiy conditions on the global stiffness matrix are applied To apply the boundary conditions 
the diagonal element of the stiffness matrix corresponding to the degree of freedom is placed as 1 and 
( onesponding value of the primary variable is substituted in the force vector Rest all the terms in the 
same row and column of that degree of freedom are set to 0 Gauss Elimination Scheme is used foi 
solving the set of simultaneous equations This gives the value of the primary variables at all the nodes 
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Post processing 

Altoi obtaining the displacementb, post processing is done to get the stiesses m each element Foi this 
f he Sttess-Sfiain Relations aie invoked and the stress tensor at each element is calculated By this slic'ss 
(('iisoi principal stiesses are calculated To find the max stress at any point the Maximum Sham Eneiqij 
Tlieoi ij (Voii-Mises Failure Criteria) is applied This gives the von-mises stress at each element Total 
(lisi)laceinent at any node is obtained by getting the vector sum of displacements m different diiectioiis 
at any node 

• MfiMinuin stiess over the clement is the maximum Von-Mises stress developed m the body 

• Ma\inuiin displacement over the element is the maximum of the total displacements occuriing a( 
(vu h node m the body 

2 8 Binary 2-D crossover 

I'Ik' inateiial domain is represented as the two dimensional domain So the crossover which operates 
(111 two diiiiciisional strings is formulated This crossover respects the shape of the parents, whcieas tlu' 
noinial onc' dimensional crossover does not guarantee to do so, hence found less effective for solving tlic' 
piohlc'in The algorithm is given m Figure 2 11 The step-wise details are given as following 

Step 1 Select two individuals Generate a random number and check whether the crossover is to be 
done' 01 not, based on the crossover probability 

Step 2 If crossover is to be done flip a com and make a decision whether to go for a, row ciossovei oi 
(olumn crossovei 

Step 3a If row crossover is to be done, for each row generate a random number and check whether 
random number is less than Ptrow/No of rows If yes, swap the rows else not 

Step 3b If column crossover is to be done, for each column generate a random number and chock 
whether random number is less than Peross/ No of columns If yes, swap the columns else not 

Thus this crossover on an average swaps at least one row or one column in the two dimensional array 
The working of this crossover will be more clear with an example 

For crossover two parents arc' chosen Let the first parent be PI and second be P2 These two have 
to be crossed with a probability Pc of 0 800 
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TIk' 10 bits paients when the shape is represented as a two 
shown li('i(' The fust low shows the column number and the first 
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dimensional an ay of size 4 x 4 aie 
column shows the low numbei The 


sK'pwise pioceduie is given lieic 


Step 1 Geneiate a landom number R Let R — 0 678 Since R < Pc, hence, the paients aie selected 
foi (lossovei 


Step 2 R = 0 185 Since 7? < 0 5 so row crossover is performed 


Step 3 P, = 

p 

JN, 

Foi 

low 

1 

p = 

Foi 

low 

2 

p = 

Foi 

low 

3 

p = 

Foi 

low 

4 

p = 


= 0 800/4 = 0 200 

0 238 Since R > Pr so this row is not swapped 
0 872 Again R> Pr so this row is not swapped 
0 156 Now R < Pr so this row is swapped 
0 333 Since R> Pr so this row is not swapped 


Step 4 Ciossovei foi this pan of paients is complete 


Aft,ei the ciossovei the childien Cl and C2 are obtained These are shown as following 
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2.9 Closure 


In this fhaptei the methodology of applying GA to the single objective shape optimization pioblems is 
(hs( uss('d 111 detail The binary strings which can be processed by the GA aie conveited into meaningful 
sh(i])es by using the connectivity mfoimation and smoothing techniques The con&tiamts aie evaluated 
1 )\ hiidmg out the stress and stiam values, by using the finite element analysis The constant stiain 
iiiangk' is used as the basic element for the finite element modelling An innovative two-dimensional 
( 1 ossovei opeiatoi is used for the ciossovei of the strings This crossover operatoi respects t lic' geometry 
and IS found very useful here To handle the constraints the selection operatoi is modified to handle 
tli(' ( onstiaints The hill climbing local search is used to ensure the global convergence of tlu' solution 
()l)t<\in(Hl by the GA For local search every time one bit mutation is carried out in a deter ministK 
in<nui('i and each time the change is accepted only if it improves the solution The combined method 
ol GA and the local search is termed as the hybrid approach 


Chapter 3 

Single objective problems and results 


In lli(' piovious chaptei, the methodology of applying GA’s to the shape design problems is discussed 
ni dc'tail and the conversion procedure of binary strings into meaningful shapes and vice-vcisa is 
( stablisluvl These shapes undeigo the finite element analysis for evaluation of stresses and strains 
In this chapter, the efficacy of the hybrid approach is demonstrated with the help of a number of single 
ol))(H five optimization problems First, all the different problems taken into consideration aie described 
riien the utility of the hybrid approach over the GA alone and the hill climbing alone will be shown 
\\ it h the help of the examples The results obtained by the use of hybrid approach on different problems 
ai(' jiK'sented and discussed The problems taken under consideration are following 

• Design of cross-sec tion 

• Design of cantilever plate with a point load 

• Design of simple supported plate for a number of different loading conditions 

• Design of hoister 

• Design of bicycle frame for static loading 

All these problems are considered for two cases 1) When the weight of the design is not considered for 
( alculation of the stresses and 2) when the weight of the design is also taken into account while finding 
stresses and strains 

3.1 Test problems 

111 this section different test problems, used to show the efficacy of the proposed approach, are described 
Foi all the test problems, a rectangular plate is taken as the basis shape for the design For the test 
pioblems following properties of the material are used 

Plate thickness 50 mm 
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Density of mateiial 7800 Kg/m^ 

\'ield strength 150 MPa 

\ouugs modulus 200 GPa 

I '(iissou s latio 0 25 

\la\ allowed Displacement 2 mm 


Foi all the test pioblems othei than the problem of design of cross-section following GA paiaineUns 
aio sc'lcH ted 


Doptdatiou Size 
t'iosso\('i Piobabihty 
Mutation Probability 
\o ol Generations 


30 
0 95 


1 /String length 
100 


To deter mine the effect of weight of the design, weight is also considered as internal loading for 
isuh problem Since with the given density, the weight of the plate is very small as compared to the 
('vtc'iiial load applied, it cannot produce much stresses Thus the effect of considering the weight will 
iiol !)(' (kvii with the normal density Hence a very large density is used, which makes the weight of 
t lie' jilatc' oi design into consideration, equivalent to the applied loading and therefore giving substantial 
( outnlmtum to the stresses and strains When weight of the structure is also considered, the density' 
IS lake'll as 1000 times the normal density, as quoted above This analysis of considering the weight 
('fiec t. ( an be useful m the applications where accuracy requirements are enormous or where the design 
IS earned out at micro level A brief description of different test problems is given as following 


3.1.1 Design of cross-section 


I'list pioblem considered is the design of the cross section which has the maximum moment of inertia 
to weight ratio The moment of inertia (MI) is calculated about its centroidal axis The distance of the 
( entioidal axis from top of the plate Yc is calculated as 


* Vi) 
* C 

Z^Wt 


(3 1) 


\\ heie 


Wt Weight of the ith element of the shape 

Distance of the center of gravity of the ith element from the top of the plate 
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This IS an unconstrained problem The weight of the element can be replaced by the area of the element 
ris tlinkness of the plate and the density of the material is constant The moment of inertia can be 
( al( ulat('(l as 

M I = Y^{w^ * [Vr - Ycf) (3 2) 

Dosiic'd latio IS given as 


Ratio = M 1 1 Weight 

I li(' plate IS divided iii the giid of 10 x 10 i e a binary stung of length 100 bits is lepiesc'iitiiig tlu' 
shape's The length and the bieadth of the plate is taken as 10 units each In this piobleni the finite 
('Ic'inent analysis is not requned as no force is applied and stress and strain values are not lequned The 
jiiobleiu IS of maximization of the MI to weight ratio for one axis 


3 12 Design of the cantilever plate 


Foi a nuinbei of engineeimg problems, structures of cantilever shape are useful , for example the hanging 
])ait of the loofs, etc The next problem taken into consideration is the design of a caiitilevei plate 
( allying an end load P = 10 kN as shown in Figure 3 1 The rectangular plate of size 60 x 100 rnm^ is 



Figure 3 1 The loading and support of the cantilever plate are shown 

divided into a grid of 60 small regular elements i e a binary string of length 60 is used to represent the 
shape The feasible shape satisfies all the constraints The optimal shape is designed for the minimum 
weight When the weight of the cantilever plate is also considered, the density of the material is taken 
siuli that the load due to the weight is comparable to that of externally applied load This requires a 
v('iy high density which may not be existing but it is used only to find the effect of internal loading 

3 13 Design of the simple supported plate 

The simple supported structures are also one of the most common structures in engineering field Few 
examples are the bases of the machines, railway tracks, bridges etc The rectangular plate of size 60 x 100 
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inur IS used foi the design This is divided in a grid of 60 elements The simply supported plate is 
(l('sigii('d loi a numbei of loading conditions as given below 

Case 1 A point load ot P = 10 kN is applied on the middle point of top side of the plate as shown in 
t h(' FiguK' ] 2 

Case 2 dist ubuted load of P = 10 kN is applied on the topside of the plate as shown m the Figure 3 3 

Case 3 A ])oiiit load P = 10 kN is applied on the middle point of the lower side of the plate as shown 
111 t li(' 1 igiiic' j 4 

Case 4 A dist ubuted load P = 10 kN is applied on the lowei side of the plate as shown in the Figure 3 5 

Foi rill t lu'se diffeient loading conditions the feasible structure satisfies all the constiamts The piobleins 
ai(' th(' iniiiiniization of the weight 

3 1.4 Design of the hoister 

Foi a uunibei of a])plications like crane hanger etc, hoister type shapes aie used for transport of mateiial 
lioiii one phu (' to another This structure has a lot of importance as by analysis of this problem 
( (' 1 1 <1111 ia( ts about the sensitivity of the shape towards increase m load can be studied This also tests 
<ip]>li( ability oi h’^ bi id approach to find shape of difficult and complicated problems of solid mechanics 
1 ic'U' some' infonnat ion about the problem like cavity information is required so that the search is guided 
lo find tlie solution having a cavity 

The rectangular plate of size 60 x 80 mm^ is used with 48 elements Thus the 48 bit long binar}' 
stung is used to represent the hoister like shape The loading and support conditions are as shown 
in the Figure 3 6 Two center elements are forced to made absent (that is, the corresponding bits in 
th(' stung aie always set to zero) The vertical load of 5 kN is applied in the middle of the cavity 
<uid horizontal loading of 2 5 kN is distributed over the length of the element For the single objective 
pioblems the objective is the minimization of the weight 

3 15 Design of the bicycle frame 

This IS a very interesting problem, which has got a practical importance, as this clearly shows the use- 
luluc'ss of this work foi the application purposes This also helps the designer to conceptualize the new 
ideas about the shapes It will be shown that this approach is helpful m finding new creative and more 
innovative solutions also Designer can also make design of his requirements by changing the constraints 
or problem definition 

The loading and supports conditions for the bicycle frame like shape is given in the Figure 3 7 In 
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i< i 2 Simple supported plate with a point 
Itiad on lop 



Figure 3 3 Simple supported plate with a dis- 
tributed load on top 


100 mm 



Figiiie 3 4 Simple supported plate with a point Figure 3 5 Simple supported plate with a dis- 

load on bottom tributed load on bottom 


this pioblem the part of the plate is made void to accommodate the front wheel A vertical load is 
applied on point A where seat is to be placed The frame is supported at two places B and C Point 
B maiks the position of the axle of the rear wheel Another point C serves the purpose of the handle 
support The filled position is marked as the location of placing the pedal assembly The shape is 
modelled as a simple supported plate and the design is conceived For this design the material yield 
si 1 ('ss is 140 MPa, Young’s modulus is 80 GPa and Poisson’s ratio is 0 25 The maximum allowed 
displacement is 5mm and the thickness of the plate is taken ^ 20 mm 
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60 mm 


Figuie 3 6 The loading and support of the hoister plate is shown 



Figure 3 7 The loading and supports for the bicycle fiame is shown 

3.2 Comparison of the results of random search with hill climbing 
and hybrid approach 

Befoie using the hybrid appioach, it is imperative to show the efficiency of the appioach considered 
foi solving the problems and advantage of the suggested approach over other potential methods In 
1 his section, this is shown by presenting the results obtained for the test problem of design of cantilever 
plate, with an objective of minimization of the weight The comparison of the results, obtained by using 
lull climbing (HC) on a solution obtained by the random search, and the results obtained by using the 
suggested hybrid approach is presented The results of the random search and GA search are obtained 
as the intermediate results 

To find the solution by GA a population of 30 individuals is taken The GAs are allowed to run for 
25 generations Other parameters of the GA are the same as given in the Section 3 1 The best solution 
found by GA is then recorded and this solution undergoes the single bit hill climbing, a local search 
meted The solution obtained after the local seardi is optimized for a single bit i e this solution can 
not be improved by the random change of any bit 
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For the hill climbing solution, first the best solution of 750 randomly cieated individuals is found 
ll('i(' 750 individuals aie taken to provide fair opportunity to the hill climbing when compared with GA 
Sni((' (F\ simulation imis foi 25 geneiations with a population of 30 so appioximatcly (25 x 30) 750 
indn idiials aic' piocc'ssod by GA before giving the best solution That s why the best solution of 750 
1 uidomlv (icMted individuals is taken foi the hill climbing The best solution imdeigoes the single bit 
l()(al hill (limbing soauh to leach the final solution 

Table' 3 i picscuits the' compaiison results for 10 diffeient simulation iims each starting with the' 
iandoin po])ulation (leated by the same landoin number 10 diffeient simulations aie earned out foi 
10 diff('ieiit e(iui-spa(('d laiidom numbeis The best solution of each simulation luii is lecoided and all 
{lu's(' good solutions colk'ctively foim a set The minimum weight solution of this set of best solutions 
IS th(' Ix'st solution and maximum weight solution of this set is the worst solution of the best solutions 
The iiK'aii median and the standard deviation of this set of best solutions aie calculated This weight 
s('t IS olitamed for the random search method, hill climbing on the best solution obtained by the random 
scviidi genetic algorithm and the suggested hybrid approach All the solutions presented m the table 
have the same' units It is important to mention that the problem is the minimization of the weight 



Random Search 

Hill Climbing 

Genetic Algorithm 

Hybrid Approach 

Best Sohitioii 

30 00 

26 00 

24 00 

20 75 

Wolst Solution 

36 00 

33 25 

32 50 

25 50 

Mc^dian Weight 

33 13 

29 38 

27 38 

23 75 

Mean 

32 83 

29 45 

27 73 

23 60 

Standard Deviation 

2 603 

2 681 

2 155 

1 261 


Table 3 1 Comparison of results for 10 different simulation runs 


From the solutions obtained for a number of test runs with different random numbers, few important 
( oiiclusions can be deduced 

• The best solution obtained by the hybrid approach in 10 different simulation runs is better than 
the best solution obtained by other approaches This is also observed that the hybrid approach 
outperforms other methods m all the cases, starting with the same initial random seed 

• The worst solution obtained by the hybrid approach, is better than the worst solution obtained 
by all other approaches This is even better than the best solution obtained by the random search 
and hill climbing on this best solution 

• The mean of the solutions is also the minimum for the hybrid approach 

• The standard deviation for the hybrid search method is lesser than other approaches This shows 
the lesser variance in the solutions obtained by the hybrid approach and all the solutions tend to 
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achieve the same optima 


• Th(' solutions obtained by GAs is the next best solution after the hybrid appioach solution This 
shows that application of GAs to solve the pioblem of shape optimization, is a good aiipioach 

• Th(' standard deviation of the GA solutions is high This is due to the lack of time piovided to 
G \s to evolve If GAs aie allowed to iiin longei the solutions come close to optima and then the' 
stanclaid deviation will also leduce 

It IS (leal fioin the lesults piesented m this section that GA can be a good appioach to solve' tlu' 
^hape (k'sign pioblem On top of GA use of local seaich helps the solution to reach the global optima 
(|ui( klv <111(1 the cpiahty of the solution is also improved This gives an idea of using the hybiid appioac h 
loi otlu'i shape d(^sign problems Moreover, the data piesented here clearly shows the efhcu'iicy of tlu' 
h\ bud appioach The mean value, the standard deviation and the best values for hybrid appioacii ai(' 
b('t t('i than any othei method Thus this approach can be said as the best approach of all the consideied 
iiK'thods for evolving the shapes 

3.3 Results for different test problems 

1 li<’ l('st piobleins aie desciibed in an earliei section In this section, the results obtained horn the 
siumlatioii luiib aic piesented The results obtained by the hybrid approach are shown in all th(' 
piobU'iiis and in some problems where the results obtained by GA are interesting or the results have 
sigiiihc ant difieiencc, intermediate results are also presented This local search on the GA solution cithei 
uupiovc's the solution oi at least the solution remains as it is But m no case the solution obtained 
1)\ GA IS degiaded by this approach This feature makes the hybrid approach more attractive as this 
.ipiuoach preserves the best solution 

3 3 1 Design of cross-section 

Foi this problem the population size taken is 100 GA simulations axe run for 100 generations with 
ri (lossovei probability of 0 95 and mutation probability of 0 01 before switching to the local search 
nu'thod The shape obtained at the end of the GA runs is shown m Figure 3 8 The local search on 
t his solution further improves the solution The final shape obtained by the hybrid approach is given 
111 Figuie 3 9 The moment of inertia to weight ratio of the best solution obtained by the GA is' 14 96 
\uiits This ratio improves to 15 72 units after the local search on this solution There is a significant 
iiiipiovement m the solution, as the shape is smoothed and the shape is similar to I-section By the 
knowledge of theory, it is known that I-section has the highest section modulus, and the same is found 
1)V the proposed approach Thus, the result is verified, as this follows the theory Since, m this case the 
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1 it;ui(' 5 8 GA solution foi design of cioss- 
-.i ( I ion 



Figure 3 9 Final solution aftei the' kx al soaic h 


tiiiinulation is such that there is no special advantage of placing of the vertical component exactly in 
I lu' niiddlo the pioposed algorithm puts it somewhere along the length of the plate That s the leasoii 
I lie soliifion docs not have the vertical bai exactly in the middle 

3 3 2 Design of the cantilever plate without considering its weight. 

I li(' next, problem is the design of the cantilever plate The loading and the support conditions aie the' 
s.iiiK' as shown in Figme 3 1 The solution obtained after the GA runs is shown in Figure 3 10 Aftei 
uMiig the local search, the solution improves and it is shown in Figure 3 11 Table 3 2 compares thc-^se 
two lestihs 


Parametei 

GA Solution 

Hybrid Approach 

Weight (units) 

22 75 

20 25 

(^max (Mpa) 

136 20 

134 68 

5max (mm) 

0 147 

0 149 


Table 3 2 Parameters of solutions obtained by GA and hybrid method 



Figuie 3 10 GA solution for cantilever plate Figure 3 11 Final solution after the local 

design search 
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'Flu' solution obtanu'd by tlu' GA has a small (hunk of the material appealing in tli(' last low ol the 
shape' drills IS cliiiimated by th(' hybiid appioach The solution obtained by hybiid appioach, is vc'iy 
siuiilai to the' paidbohc shape, which is the solution with unifoim stress Interestingly, this appioach 
('liiiimates the two lows fiom the solution, one from the upper side and one from the lower side, showing 
I li(' potential of finding the optimal use of the material This is clear that the hybrid appioach is able to 
uiipiove the solution Moie mteiesting observations are the values of the maximum stiess and strains 
i;iv('n 111 Table 3 2 The solution obtained after using the hybrid approach has lesser maximum stiess 
Ilian the' GA solution, but there is slight increase m the deflection This is due to the re-orientation of 
I he shape The results show that the solution obtained is near to the stress constraint This mean that 
the stress constraint is more important constraint for this problem 

3 3 3 Design of the cantilever plate when the weight of the design is also considered 

rii(' ('xtc'iiial loading conditions remain the same as previous problem Only difference is that the load 
of the structure itself is also accounted as a distributed load, foi the calculation of the load The solution 
obtained aftei the GA runs is shown m Figure 3 12 After using the local search, the solution improves 
,111(1 it IS shown 111 Figure 3 13 


Par ameter 

GA Solution 

Hybrid Appioach 

Weight (units) 

27 50 

24 25 

^max (Mpa) 

135 26 

142 22 

5mai (mm) 

0160 

0 172 


Table 3 3 Parameters of solutions obtained by GA and hybrid method 



Figure 3 12 GA solution for cantilever plate 
dc'sign (when weight is also considered) 



Figure 3 13 Final solution after the local search 
for the cantilever plate(weight is considered) 


Hcie the shape obtained by GA, has a crossed rib This rib can reduce the stresses But moie 
ml ('icstmgly this iib is eliminated by the local search in order to save material This causes a slight 
nuK-^ase m the stress but the solution remains feasible Though the solution obtained by the hybrid 



ai)pioadi IS veiy similar to the parabolic shape, which is the solution with uniform stress, yet the 
solution is significantly different from that when the weight of the plate is not considered The solution 
(ihi.uiH'd by the hybiid approach make better use of material as shown m Table 3 3 

1 3 4 Design of the simple supported plate when a point load is applied on the top 

\('\r ]>i()hl('in taken into consideration is the design of the simple supported plate, when a point load 
IS ai)])h('d on the topside of the plate m the middle as shown in Figure 3 2 

Case 1 Weight of the designed plate is not considered 
1 lu' solution obtained by the hybrid appioach is shown in the Figure 3 14 This solution is sinnlai to 
I he solution obtained after GA luns The numeiual data for GA solution and the hybiid appioadi is 
'a,i\('n 111 Table 3 4 


Parameter 

GA Solution 

Hybrid Approach 

W eight (units) 

20 25 

19 50 

^max (Mpa) 

49 43 

50 57 

Smax (mm) 

0 021 

0 022 


lalile 5 4 Paiaineters of solutions for simple supported plate with point load on the topside (Case 1) 




Figinc' 3 14 Solution obtained after the hybrid Figure 3 15 Solution after the hybrid appioach 

ripiiioac h foi the simple supported plate when a for the simple supported plate with a point load 

lionit load is applied on topside at the top (weight is considered) 

The difference in the weight of the solution obtained by the GA is not much reduced by the solution 
obtained by the local search Interestingly, the solution obtained by the hybrid approach eliminates all 
I li(' mateiial other than the arms joining the loaded point to supports The stress and strain values aie 
v('iy small This implies that the thickness of the plate can be reduced which will further reduce the 
weight and make more effective use of the material It has been observed that for a very small range of 
htness there are many solutions which have very httle difference m the fitnesses but the solution looks 
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different from each other This is due to the reason that the search space is discreet, with continuous 
sc'aich space this problem is less likely to happen The shape obtained somehow depends on the initial 
laiuloiii population Eventually, if the simulations are allowed to run longei, it may happen that due 
to genetic drift the solutions may go to one solution But this may require veiy-veiy long time befoie 
t his phenomena is observed 

Case 2 Design of the plate when the weight of the design is also considered 
I li(' jilate eonsidcied in the pievious section is redesigned for the case when the weight of the plate also 
( ontiibutes to the loading in addition to the externally applied load The shape obtained by the hybiicl 
appioach is shown in Figure 3 15 


Parameter 

GA Solution 

Hybrid Appioach 

Weight (units) 

20 00 

19 50 

(^n}ax (Mpa) 

89 44 

87 42 

Smai (mm) 

0 048 

0 048 


Table 3 5 Parameters of solutions foi the simple supported plate when a point load is applied on the 
topside and weight is also considered 

This solution obtained by the hybrid approach is slightly different in weight and the configuration from 
the solution, when the weight is not considered The solution has same weight and moie stresses than 
the case when the weight of the solution is not considered, as the loading is increased now The weight 
ol the solutioTi can be reduced by using the thiimer plate for the design The result obtained by GA is 
very similar to that by hybrid approach as shown by the numerical values from the Table 3 5 

3.3.5 Design of the simple supported plate when a distributed load is applied on 
the top 

111 this problem, the distributed load is applied on the top of the plate The loading and the support 
( oiiditions are shown in Figure 3 3 To make the solution geometrically feasible, all the top row elements 
must be present 

Case !• The weight of the plate is not considered 
The solution obtained by the hybrid approach for the case when the load of the plate is not considered 
lb shown in the Figure 3 16 

The numerical data presented in the Table 3 6 shows that the solution obtained by GA and by the 
hybrid approach is same The interesting thing about the solution is thg-t the left arm is straight joining 
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Figuu' 3 15 Solution alter the hybrid approach 
when the distiibuted load is applied on the top 
()l simple suppoited plate 


Figure 3 17 Solution foi the design of the sim- 
ple supported plate with distributed load with 
its weight also consideied 


Par ameter 

GA Solution 

Hybrid Appioach 

Weight (units) 

23 75 

23 75 

CTmax (Mpa) 

32 91 

32 91 

6max (mm) 

0 016 

0 016 


Table 3 5 Paiaineteis of solutions for simple supported plate with distributed loading when weight is 
not consideied 

the lop low with the support node The right arm is slanted This is also clear from the numerical data 
ol the (h'sigii that the thickness can be reduced as the stress and strains for this design are too less than 
I li(' <illow('d value The reduction in the thickness helps in the optimal use of the material The loading 
( an also be increased to make optimal use of the material 

MoK'ovei, the design does not seem to be optimal as the loading is symmetric, hence the solution 
obt, allied must be symmetric This can happen due to the msulBcient time available for the evolution 
thiough GA for this problem The solution undergoes local search with one bit change at one time This 
may not be sufficient to get the global optimal solution m this particular case If the solution obtained 
by the GA is tried to be optimized with 2-bit flip or more bit flip simultaneously, better solution can 
be obtained, but the computational effort required is more 

Case2. When weight of the designed plate is also considered 
The loading conditions now include the weight of the design, as the distributed load The final solution 
obtained is shown in the Figure 3 17 

The numerical data of the solution as presented in the Table 3 5 This clearly shows there is no 
improvement in the solution obtained from the GA solution Here the solution is having its right arm 
straight And the left axm is somewhat slanted The slant angle is more than the case when the weight 
of the solution is not considered The same results are obtained for the many simulation runs with 
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Parameter 

GA Solution 

Hybrid Approach 

Weight (units) 

24 25 

24 25 

^max (Mpa) 

106 48 

106 48 

6mai (mm) 

0 087 

0 087 


J 7 Paianieteis of solutions obtained foi the design of simple supported plate with a distiibuted 
load applied on the topside of plate when weight is also considered 

(liftc'ienl landoni numbeis The solution requires more GA runs to reach the global optima 

3 3 6 Design of the simple supported plate when a point load is applied on the 
bottom-side 

I’lie simple supported plate is designed foi the case when a point load is applied on the bottom side as 
shown in the Figuie 3 4 

Case 1 When weight of the plate is not considered 



Figuie 3 18 Solution after the hybrid approach Figure 3 19 Solution after the hybrid approach 

loi the simple supported plate with a point load for the point load on bottom side when weight 

on the bottom side is considered 


The solution obtained by the application of the hybrid approach is shown in the Figure 3 18 The 
numerical data for the problem is given m the Table 3 8 


Parameter 

GA Solution 

Hybrid Approach 

Weight (units) 

10 00 

10 00 

CTmax (Mpa) 

136 67 

136 67 

^max (mm) 

0 081 

0 081 


Table 3 8 Parameters of solutions obtained for the design of simple supported plate when a point load 
IS applied on the bottom side of the plate 


The results shown in the Table 3 8 clearly show that the GA finds the global optimal solution This so- 
lution IS the global optimal solution as the shape caimot remain feasible with lesser number of elements 
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'I’lns solution has mteiestmgly eliminated all the four rows from the top without having any pioblom 
luloiination To connect the node where the load is applied to the support nodes, the path taken is 
the shoitcst path So the solution becomes the global best solution Thus the hnal solution obtamc'd 
iluougli the hybiid appioach is the same as the GA solution 

Case 2 When the weight of the designed plate is also considered 

1 li(' same' plate is ledesigned foi the case when the weight of the designed plate is also (oiisideied The 
K'siilt obtanu'd by the hybrid approach is shown m the Figuie 3 19 The numeiical data foi the piobk'iii 
IS piesented in the Table 3 9 


Pen ometer 

GA Solution 

Hybrid Appioach 

Weight (units) 

13 75 

13 75 

O'nioi (Mpa) 

114 56 

114 56 

5nia, (mm) 

0 054 

0 054 


Table 3 9 Paiametcis of solutions obtained for the loading on the bottom side when the weight of 
simple supported designed plate is also considered 


H('ic' t he shape' is consideiably different from the result m the previous case As the loading is increased 
diu' to the weight of the design, the previous shape does not remain sufficient to sustain the load So 
ilu' shape is (hanged More interestingly, the solution has got two cavities Here also the result does 
not show any iinpioveiiient by the local search It means GA is able to find the optimal solution in 
tins (as(' The impioveinent is only possible by reduction in the thickness Table 3 9 gives the value of 
Stresses and stiaiiis for the optimal solution obtained by GA and by hybrid approach 

3 3 7 Design of the simple supported plate when a distributed load is applied on 
the bottom 

The loading and the support conditions are shown in the Figure 3 5 The problem requires all the 
elements of the last row to be present in the solution to make it geometrically feasible 

Case 1 When weight of the plate itself is not considered 
I'he design is earned under the action of external loads only The solution obtained by the hybrid 
appioach foi this problem is shown in the Figure 3 20 

The numerical data for the solution is given in the Table 3 10 

The shape obtained by the hybrid approach eliminates the top rows to reduce the weight Since all 
the elements in the bottom row must be present, to make the geometry feasible and there is no extra 
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1 J 20 Solution hs the hybiid appioach 

w luMi the' (hstiihuted load ls applied on the bot- 
loin si(l(' of simple supported plate 



Figure 3 21 Solution after the h\biid approach 
for the srmple supported plate with the dis- 
tributed load on bottom side when wc'ight is 
considered 


Pa7 ometer 

GA Solution 

Hybrid Appioach 

Weight (units) 

10 00 

10 00 

^niax (Mpd) 

78 20 

78 20 

^ma i (min) 

0 052 

0 052 


lablc' 3 10 Parameters of solutions obtained for the design of simple supported plate with a distiibued 
loading applied on the bottom side when the weight is not considered 

('Ic'inent m the solution, this solution can be said a global optimal solution The data shows that GA 
IS able' to find t.he global optimal solution Since no additional element is present there is nothing to 
impiove by the local search, hence the same solution is retained by the local seaich The weight can 
only be reduced by reducing the thickness 


Case 2 Weight of the designed plate is also considered 
The solution for this case is shown in the Figure 3 21 The solution is very interesting as this shows the 
tiond of placing the material for this problem when the loading is increased The numerical data for 
the' solution show that the global optimal solution for this has more weight than the solution obtained 
loi tlie case when weight of the designed plate is not considered 


Parameter 

GA Solution 

Hybrid Approach 

Weight (units) 

14 00 

14 00 

^max (Mpa) 

142 69 

142 69 

8max (mm) 

0 069 

0 069 


Table 3 11 Parameters of solutions obtained for the design of simple supported plate with a distributed 
loading applied on the bottom side when the weight is also considered 


Since increased loading is present, hence the solution requires more number of elements to sustain the 
load The result obtained by GA does not improve by the local search 
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3 3 8 Design of the bolster 

1 li(' loaclmg c onciitions and the supports are shown in the Figure 3 6 This problem is a veiy inteiesting 
pinhlc'ui as the ('\p<'( tod solution have a cavity The simulation results after completing the piocess oi 
C \ and th(' lo( al so.iuh aio given m as following 

Case 1 Weight of the hoister itself is not considered 
th(' v('ight of the designed hoistei is not considered The solution obtained is shown m th. 

I igiiu' 1 22 



1 iguK' 3 22 Solution altci the hybrid approach Figure 3 23 Solution after the hybiid appioach 

loi th(' di'sign of hoistc'i for the design of the hoister when weight is con- 

sidered 

Th(' solution IS veiy inteiesting m the nature, as a perfect hook is obtained The solution is not a 
( losc'd loop solution, which has lesser stresses but the weight is more The result eliminates the leftmost 
.ukI iightmost (oliniiii and the bottom most row This hoister can be used for carrying the goods as in 
the ciaiics etc This solution is the least weight solution Most interesting observation for this pioblem 
IS that foi the same fitness, three different solutions were obtained This means the configuration of the 
solutions aie dependent on the initial random population The numerical data for the design consideied 
IS given ill the Table 3 12 


Parameter 

GA Solution 

Hybrid Approach 

Weight (units) 

12 50 

11 00 

CTmax (Mpa) 

67 74 

68 77 

^max (mm) 

0 072 

i 

0 081 


Table 3 12 Parameters of solutions obtamed for the design of the hoister 


The lesult in the Table 3 12 shows that the solution obtained after local search is improved 
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Case 2 When weight of hoister is also considered 

\o\\, the same hoistei is designed when the load of the designed hoister is also consideied as having 
s.iiiK' t on tlu' loading The lesult obtained is shown in the Figuie 3 23 Now the lesult is not hkc' 
I hook Instc'ad a small difloience in the conhgination is obseived 


Parameter 

GA Solution 

Hybiid Approach 

Weight (units) 

15 75 

11 00 

(^max (Mpa) 

76 43 

88 82 

5,noL (mm) 

0 048 

0 087 


I 1 1 5 Paianiet eis of solutions obtained foi the design of hoistei when the weight is also c oiisideu'd 

llu' Table 3 13 shows the numeiical data obtained foi the design problem The coiupanson of (Ik' 
G\ K'sult and the final solution data shows that the solution obtained aftei the local seauh is not 
oiih inipiovc'd in the weight , but also it leairanges the material such as the maximum stiess and tlu' 
luaMinum deflection aie also reduced Most interesting thing about this solution is that this suhitiun 
also has got the weight same as the solution m the case when the weight of the designed plate is not 
(onsideied But the configuration has now changed This geometry has got the advantage of having 
lessen moment as compared to the previous case solution This shape has got both the aims caiivmg the' 
lion/oiilal load are connected to the body so these tend to neutralize the effect of each other inaking 
I 1k' shape' more stable This change in configuration helps m reducing the stresses when the weight of 
I 111 ' (U'sigii IS <ilso accounted for the loading 

3 3 9 Design of the bicycle frame 

Tins IS the most interesting problem as this is closest to the real world problem The solution of this 
pioblem shows the ability of the suggested approach in solving the real world problems The support 
.lud t,hc' loading conditions are shown in the Figure 3 7 

Case 1 Weight of the frame is not considered 
Fust, the design exercise is carried out when the weight of the design does not account for the loading 
The best solution obtained by the hybrid approach is shown m the Figure 3 24 


Parameter 

GA Solution 

Hybrid Approach 

Weight (units) 

35 25 

28 50 

^max (Mpa) 

38 03 

32 13 

hraax (mm) 

0 157 

0 118 


Table 3 14 Parameters of solutions obtained for the bicycle frame design problem 
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I ionic i 24 Solution dftei the hybud appioach Figuie 3 25 Solution dftei the hybnd appioach 

Ini till' design of bicycle flame foi the design of the bicycle fiaiue when its 

weight IS consideied 

I 111' iK'st lesult IS veiy much different hoin the conventional shape of the frames This has got a stiiu- 
niic which has eliminated the aim connecting the paddle support block and the handle support nodes 
1 li(> shape of the bicycle with this frame is shown in the Figure 3 26 This solution can be a potential 
( andidate of design in the search of the new models of the bicycle But this designed is required to 
1)(' (lic'ckcd for dynamic loading conditions The numerical data from the Table 3 14 shows that the 
St losses of the solution are very small Thus the material can be taken from the lateral diiec tioii without 
alh'ctmg the shape The data also shows the reduction m the weight of the solution obtained by the 
GA l)y the local search 

Case 2 When weight of the frame is also considered 
Hc'ic' the design is nearer to the practical conditions The weight of the solution is also considered, but 
the' density of the material used is increased by 200 times than that taken in the case when weight of 
rlic' design is not considered This is done to have the internal load due to the weight of the same order 
ol the external load 


Parameter 

GA Solution 

Hybrid Approach 

Wetght (units) 

42 00 

31 25 

CTmax (Mpa) 

86 39 

79 53 

^max (mm) 

0 709 

0 910 


Table 3 15 Parameters of solutions obtained for design of bicycle frame when weight is also taken into 
( oiisideration 


The solution obtained by the hybrid approach is shown m the Figure 3 25 The solution has more 
weight than the previous case Table 3 15 presents the numerical data for the design The stress and 
slaam values are more for this case than the previous case but still the solution can be improved by 
the reduction in the thickness It is very clear from the numerical data obtained that the local search 
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^Imimatos some pait of the material which does 
li(> (lottc'ction of the design has inci eased 



1 in UK' 2b Bicycle shape foi the fiame when 
li(‘ wc'ight IS not coiisideied 


participate in the stress sharing In the process 



Figure 3 27 Bicycle shape toi the haine when 
the weight is consideied 


not 


3 4 Closure 

1 Ins (ha}:>tei deals with the description of the various shape design problems solved to show the efficacy 
[)l the appioach Then the results obtained by applying the hybrid approach to solve the tesi pioblcms 
,u(' ])ieseiited and discussed The comparison of the local search method with the hybiid method 
IS iisc'd to show the woith of the pioposed approach The results support the concept of using the 
li\l)iid appioacdi to solve the shape design problems The results obtained for the single objective 
()j)i imi/at 1011 levc^al some inter esting facts The stress constraint is found to be the active coustiamt in 
all the' pioblems The solutions obtained by the considered approach do not have any control on the 
( oiifigiuation It is possible to obtain different shapes which possess the same fitness value Thus the 
1 ('suits provide sufficient motivation to solve the problems as the multiple objective problem Then the 
vaiiation in the weight as well as the effect of configuration on the design can be studied more easily 
1 he' tunid of weight change and optimal method of placing the weight in the design, such as it is used in 
th(' most efficient way, can also be obtained by the study of the problems for the multi-objective caises 
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Chapter 4 

Multi objective GA for shape design 
problems 


111 lh(' s(H 011(1 cliaptei, it is discussed how GA’s can be applied for a single objective piolilem In tins 
( liri})t ('1 It V' ill be shown how the Pareto-optimal solutions can be obtained for the shape design pioblcniis 
with iiioK^ than one objective functions This study deals with the case when there aie two ob]e(ti\('s 
ol ( oiiHu ting iiatiue The pioblem differs from the single objective problem in the sense that the uiun- 
l)('i of ol)]e( fives is moie than one, and these are of conflicting nature Then no single optimal solution 
( an l){' said t he best solution, instead a number of solutions which he on Pareto -optimal Front aic found 

h\n thc' classical methods there is a single way of solving these types of problems— convert tho 
mult i-()h](K tive problem into a single objective problem There are a number of ways of doing this 
wcnghtc'd sum appioach, ^perturbation, mm-max method, goal programming and others But none of 
thc's(' inch hods aic^ reliable as the solution largely depends on the parameter settings to formulate the 
single objective problem Besides this the classical methods give one solution at a time, so to get a 
nuinlxn of solutions, many simulation runs are required Multi-objective GA’s do not face these piob- 
h'lns while dealing with a number of objectives as well as a population of individuals A numbei of 
nmlti-ohjoctive optimization techniques [6, 15, 21, 24, 29, 32, 38] are suggested and it has been shown 
that these evolutionary studies has the potential to achieve the Pareto-optimal front as well as they can 
leach the global Pareto-optimal front [37] 

In this chapter a hybrid method based on the combination of evolutionary algorithms and the local 
s(vu( h method is suggested to find the optimal shape design problems The number of solutions obtained 
attei the evolutionary search is more than required, so these are reduced to a practically suitable numbei 
by using a clustering method 
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4.1 Multi-objective optimization and pareto optimality 


111 (\ ninlti-ob] 0 ( tivc optiimzdtion thcic fiio nioie than one objectives If these objectives have the opti- 
mal solutions (iistaiit to ca<h othci the objectives aie said to be conflictmg in natuio Multi-( iiteiion 
()[){ iini/at ion with su< h (onfiicting objectives lesults in a set of optimal solutions Many solutions aic' 
opinnal because no solution can be considered bettei than othei with lespect to all the objectives This 
solutions set is bc^ttcr known as Paieto- optimal solutions 


Definition 2 Pa) vto-optnval Fiont is the set of bolufions which are equally good as coaipaied to aaij 
i)fli( I sohiiioi) of the set ? e ariij solution on the front can not be selected when conipaied with oihet 
iiitlioiii b( niq biased towards other MathematLcallg %t can be given as 
1 (le( ision vector i E Xf rs said to be non-domrnated regarded to a set A C Xj if 

e A ay X 

// if IS (lea) within the context which set A is meant, it is simply left out Moreouei, xis said to he 
PuK to-opfnnal ifxis non dominated regarding Xj 

iM.t’UH' 1 1 givers a bcttei illustration of the concept of Pareto- optimality through a example pioblem of 
tiadc^-ofi between two conflictmg objectives namely, cost and the accident rate, both to be minimized 
The point A lepiescnts the solution with a very small cost but highly accident prone Another solution 
b u'pu'seuts the solution very safe but very costly Here both the solutions are good with lespect to 
tlu' one objective, but worse m other objective There exist a lot of solutions (like solution D) which 
also belong to the Pareto-optimal set and have equal importance to a unbiased decision maker All 
llicase solutions (lying on the thick solid line) are known as the members of Pareto- optimal Front and 
this solid line is called Pareto- optimal Front 

Accident rate 


Pareto-Optimal Front 

Figure 4 1 The concept of pareto-optimal solutions is illustrated 




45 


i 



4.2 Dominance 


Tins IS a \nv nnpoita.ut concept so far as the discussion on multi-objective problems is conceined This 
loiiiis th(' l)<u k 1)0110 of the ranking strategy For constrained and unconstrained dominance dohnitioii 
\,ui('s Both o{ the clofiiiitions aie given as following 

4 2 1 Dominance for unconstrained problems 

Dehnitioii 3 If the n cue N objective functions and two individuals, “a’' and “b” are compaied then 

Cl >- 6 (a dominates b) i//(a) > f{b) 

f{a) > f{b) means that \/ieN Ua) > Mb) 

01 

a b (a. weakly dominates b) if/(a) > f{b) 
fio-) > f{b) means that 
Vt e iv Mo) > Mb) 
and Ai £N fM) > Mb) 
or 

a ^ b (a IS indifferent to b) if f{a) ^ /(6) 

A f{b) t f{a) 

In suuiilo wolds i1 c an ho said, that “a dominates b” if and only if for all function values, the function 
vahu' ol a is oithoi better or at least equal to that of function value of 6, and at least one function value 
of 0 IS better than b If all function values of a are equal to those of b or if in some function values a is 
l)ot,toi and for some function values b is better, a and b are said to be indifferent to each other 
This ( an bo uiidoistood with the help of the Figure 41 It is observed that the solutions which are not 
on the Paioto-optimal front (like solution C) also exists If it is compared with solution A, according 
to defiiiitioii 3 both are indifferent But if it is compaied with another solution D, it is found to be 
dominated by solution D Since for a solution C, there exists a solution D in the search space, which 
dominates it, this solution cann’t be m Pareto-optimal Front 

4 2 2 Dominance for constrained problems 

In the multi-objective cases where the constraints are present, definition of “dominance” is redefined to 
accommodate the constraints 

Definition 4 The definition of dominance when more then one constraints are present is given according 
1o the following conditions 

Condition 1: If individual A is feasible and individual B is not then A dominates B 
Condition 2 If individual B is feasible and individual A is not then B dominates A 
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Condition 3‘ If both of the individuals are infeasible then the individual with lesser violation of con- 
^fiawt(h) dominates the other 

Condition 4 If both of the individuals are feasible or equally infeasible, dominance is decided acjoidmij 
to 1h( definition I 

4 3 Representation scheme 

1m)i limit i-ol))e( five pioblems the lepiesentation scheme (Section[ 2 2]), finding connectivity (Se(tion[ 2 2 2]] 
smoothing opoiation (Section[ 2 2 3]), material selection (Section[ 2 3 3]), and finite element anahsis 
(S('( tion[ 2 7]) rue the same as discussed m second chapter Here the fitness is represented as a \c( toi 
as multiple ol))ec tives aie present instead of a single objective So the selection method is also diffeieiit 
lo sohc' siuli pioblems we deviate from simple GA and the specialized algorithms— inulti-objei tni' 
gcmetic algoiit hms (MOGA) — which can handle multiple objectives aie used Foi this purpose a specifn 
MOGA Non -dominated Sorting Genetic Algorithm -II (NSGA-II) [6] is used 

4 4 Multi- objective problems 

Foi multi-ob]e( tivc pioblem of shape design problem, two minimization type objectives are consideiecl 
riiKH' (()iisU<mit need to be satisfied to have a safe design 

Objective functions 

• Fust objective' is weight of the object It is defined as 

Weight = Area * Thickness * Density (4 1) 


This IS to be minimized 


• Second objective is stiffness It is defined as 

Stiffness = Y ? — r 

[^maxtmuTnl 

This objective is to be maximized This can also be formulated as minimization of maximum 
displacement Then both objectives are minimization type 

Both of these objectives are of conflicting nature, as more material means more weight but loosely 
saying moie stiffness and lesser material means lesser weight but low stiffness So there is no single 
optimum solution, instead it gives a set of optimum solutions which he on the Pareto-optimal front 
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Constraints 


SuiulfU to tho singlo ob]ective case three constraints aie theie but there is a little diffoience lu the 

( oust I amt violation (oiulitions 

• Gc'oiiu'tn Coustidiiit This IS same as desciibed in Section 2 3 2 

• St loss Coiistiaint This is also same as described m Section 2 3 2 

• DispLuoiiu'iit Coiistiaint Here we put a uppei limit on the allowed displacement If the' nia\ 
(lisplac c'liK'iit IS inoie than the specified value the solution is not accepted But now tins constiaint 
does not (ontiibutc m finding the enoi 

4 5 Non-dominated sorting genetic algorithm-II 

This algoi ithiii has the features required for a good MOG A It is shown elsewhere [7] that it can niaiiit am 
t lie (livc'i sity on the Paieto-optimal front as well as can converge to global Pareto-optimal Fiont Besides 
this the time order complexity of this algorithm is shown as 0{MN^) [7] The steps involved m the' 
algoi it Inn aie gi\en here 

Step 1 Initialize' the population 

Step 2 Calc ulate the fitness 

Step 3 Rank the population using the dominance definition defined m Section 4 2 or m Section 4 2 2 
as applicable for non-dominated sorting [32] 

Step 4 Calculate the crowding distance [Section 4 5 1] 

Step 5 Do selection using crowding comparison operator [Section 4 5 2] 

Step 6 Do crossover, mutation and generate the intermediate children population 

Step 7 Combine the intermediate population and parent population and do non-dommated sorting 
and find crowding distance 

Step 8 New parent population is formed by accepting solutions from the first (best) non-dommated 
front and continuing to other fronts successively till the new population exceeds the number of 
individuals as the parent population had 

Step 9 The solutions from the last accepted front are sorted according to the crowding distance and 
as many individuals are selected, which make the population of new parent population same as 
old one 
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step 10 Repeat Step 5 
4 5 1 Crowding distance 

1 11 okU'i to piescive the diveisity along the fiont, the selection of individuals which aic vciy neai to eadi 
otlicu must be discoiuaged To do the same this concept of crowding distance is used The ciowdmg 
(list aiK e IS nieasuied as the distance of the biggest cuboid containing the two neighboiuing solutions oi 
ih(' scUiK' Pau'to-optiinal fiont in the objective space as depicted in Figure 4 2 The ciowding distance 



Figiiie 4 2 The crowding distance calculation is shown 

IS (akulated by fiist soitmg the solution set with respect to each objective function and finding the 
(hstaiicc between the immediate neighbours The solutions on the ends aie assigned a large value of 
(lo’wdmg distance This helps m preserving the extremities of the solution set Then the sum of these 
(hst-ances foi each solution is called as crowding distance of that particular solution For more details 
about the ciowdmg distance lefei to [6] 

4.5.2 Crowding comparison operator 

Since the diversity among the solutions is important, a crowded comparison operator is defined using a 
icTation -<n as following 

Definition 5 Solution % is better than solution j in relation -<n if (irank < Jrank) or ((irank = Jrank) 
mid (idhstance ^ Jdistance)) 

wheie i distance IS the crowding distance of the ith solution and irank is the rank of the ith individual 

4-5 3 Selection 

In the single objective problems the selection is based on the basis of the fitness of the individuals but m 
multi-objective optimization problems there are more than one objective functions, so selection cannot 
bc"* done on the basis of fitness as both solutions may be better than each other m different objectives 
Selection in multi-objective problems is done according to crowding selection operator described m the 
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pK'vious section Better rank individuals are given preference over the individual with worse rank of 
t he two If both individuals have the same rank, the individual with more crowding distance is selected 

4 5 4 Elitism 

Chtisin IS a way of pieseiving the good solution found so far This helps the evolution strategy to 
(oiivc'igc' fastei to the global Paieto-optimal front To employ the elitism, old and the population 
'^('lu'iated aftei selection ciossovei and mutation axe combined This combined population is laiiked 
uid tins population is halved, selecting the best individuals This way all the good individuals au’ 
pu'sc'ived 

4 6 Hybrid methods 

riu' evolution stiategy selected for solving multi-objective problem is found to peifoim veiy well on 
till' tc’st pioblems, which are of very complex nature In this section, a hybrid stiategy combining the' 
c'M'lutionaiy method with classical method, is suggested to make the appioach more piactical bv 

1 ensuimg convergence to the Pareto-optimal front 

2 ' educing the size of the obtained non-dommated set 
4 6 1 Better convergence 

111 a leal world problem the knowledge of global Pareto-optimal front is usually not known Though 
NSGA-II has shown the ability to find the global Pareto-optimal front for a number of test problems 
[7 8] but it IS always helpful to enhance the probabihty to reach the true global Pareto-optimal front by 
losing the a local search method at the rear end of the algorithm The local search algorithm requires 
single objective function, hence a weighted sum of objectives is used to convert the multiple objectives 
into a single objective The single objective is defined as 

M 

(4 2) 

where 

/j(x) objective function 

wf weight corresponding to the objective function 

In the equation 4 2, if the scaled function values are not used and if the magnitude of the function 

values are not of same order, the obtained sum may be biased towards one function though it is not 

desiied and intended To avoid this difficulty the fitness function values are scaled between 0 1 and 1 0 
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S( almg of the function values is done according to the following formula 


T scaled _ q , if^ 

^ ^jmax __ 

w li('l ( 

jsiahd Scaled value of ith objective function 
f I Value of zth objective function 

idiniinum value of ?i/, objective function over the population 
JcIaMinum value of in, objective function over the population 


(4 3) 


4’h(' Inggc'st obstacle to this single objective pioblem is m finding the weight vectoi as it is obvious 
Ik nil flu' ('(juafiuii 4 2 that the weights affect the sum of fitness Any setting of the weight vectoi will 
livul to a (hfic'ient optima There can be a number of ways of finding weight vectoi like giving equal 
iiiijioitauc e to each function value or let the user specify the weights etc But all these methods siiffci 
hoiii tlic' vc'iy fundamental problem of requirement of the problem information before hand Ifthe weight 
IS not chosc'ii properly, it may mislead the search and guide to any local optimal point The pioposcd 
mc'thocl of weight calculation is elaborated in Section 4 6 4 Since this method of weight calculation 
givc's the' weights different from the user specified weights so t hese are called as pseudo- weight vector 


Ojitmii/af 1011 of the weighted sum of the fitnesses, F{x) is done using the hill climbing local search 
iiu'f hod [Sec tioii 4 6 3] This optimization will produce a Pareto-optimal or a near Pareto-optimal solu- 
tion as illusf, rated lu Figure 4 4 This is true for the convex Pareto-optimal regions Since independent 
loca’ search inc'thods are tried for all different solutions of the Pareto-optimal set, it may happen that 
all the' solutions obtained by local search are not non-dominated to each other Thus, once again a 
uoii-dominated sorting is performed to obtain the non-dommated set 


4 6 2 Reducing the size of non-dominated set 

hi the multi-objective problems, a number of solutions are obtained, which are represented as different 
jioiuts on the Paieto-optimal front These gives a number of solutions to the designers or the designer 
will have' a lot of choice, but practically it is not feasible for the designer to take a look on all the possible 
solutions So the number of solutions must be reduced to a number on which, a designer can concentrate 

10 decide a final solution according to his/ her requirement It may be a question here — When so many 
solutions are not wanted, why to obtain so many by GA, why not reduce the population size‘s The 
answer to this is that it can not be done GA’s work with the population To find the global optimum 

11 lequnes a optimum level of population [17] 

Very small population will not have the potential to search ahead and reach the global optimum or 
it may prolong the convergence It will converge very early to a local optima K very large population 

,-,..133625 



IS ( oiisicleied oui search will reach the global optimum but it will be consuming a lot of time as this 
pjo(('ss ks computationally very expensive due to finite element analysis of all the shapes and then 
ili(' nou-doiuiiiated soiling Besides this again a lot of solutions will be pieseiit, which pose the same' 
piobknu to the de&ignei of having too many solutions 

To i(h1u((' the numbei of the solutions, some solutions, which can be improved oi which foims a 
( liisU'i ( an gi\(' way to othei solution which aie good as well as have diveisity in the difteient objectives 
1 Ins ('\('i(ise IS (allied out to some c'xtent by the strategy of local search and then mainly clusteimg 
,S('( lion 4 0 5] takers caie of it 

The' complete pioceduie of the pioposed hybrid method is shown m Figuie 4 3 Staitmg fiom the 
piohknn the evolutionaiy stiategy yield the primary solution set This is piocessed by the local seaich 
iiK'thod winch leads to a bettei converged set of solutions Non-dominated solutions aie extiacted fiom 
tins scy Clust('img leduces the number of the solution to a number required by the designei 



Figuie 4 3 The proposed hybrid procedure of using a local search technique, a non-dommation check 
and a clustering technique is illustrated 

4 6 3 Hill climbing 

The Hill climbing method for the multi-objective process is given here This is a method of finding the 
optima starting from a initial guess This initial guess is provided by the solution obtained by MOGA 
urns This IS more like a steepest descent search, trying to find the solution in the global basin, when 
it IS expected to be near the global optimum The proposed method is as following 

1 Select an individual and calculate the unrefined weighted sum of the scaled fitness of this solution 

2 Mutate one bit each time m the individual string 
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3 Extidct the shape from the new string 

4 Cakulate the stiesses, displacement and fitness values for this new solution (shape) 

C',il( ulatt' the wcnghtcd sum of the scaled fitness 

() Coinpaie the uiiiefined sum and the new obtained solution If theie is impiovenient iii the now 
solution accept the change else reject it 

7 111 cas(> ol icjcH tioii lestoie the pievious bit value 

8 Rc'pc'at tlie piocess till all the bits aie mutated once 

9 Calc ulatc' the total number of bit changes If theie are no bit changes, teimmate the hill climbing 
piocess foi this individual, else repeat Step 2 , from the first bit of the individual stung 

10 11 all the individuals have undergone this hill climbing process, terminate the hill climbing cdsc' 
ic'pc'at Step 1 

The' Figuic 4 4 illustiates the basic phenomenon of hill climbing The basic idea behind this cxeicisc' 
IS that once the solutions close to the true Pareto-optimal front are obtained from evolution stiatcgy 
tlu' local sc'aich method with different emphasis of objective functions can be used as an attempt to 
icMch the tiiic' Paieto-optimal front 

Tins IS also of the importance from the point of view of reduction in computational effoit It may 
hapiic'ii that the' seaich algorithm is near the true Pareto-optimal front but it may requiie inoie nnmbei 
of gc'nc'idf.ion to leach the global front, making the use of evolutionary method alone computationally 
c'xpensive Heic we can stop the process of evolution at this step and let the local search method take' 
it ovci to find the global Pareto-optimal front speedily 

4 6 4 Weight calculation 

Heic', a stiategy similar to that proposed by Deb, 2000 [10], which eliminates the requirement of weight 
deteimniation from the user, is used In fact, this methodology determines the weight by using the 
fitness values of the solutions, hence ehmmating the problem of weight determination by the user 
^Vclghts are calculated on the basis of their function values Two methods of weight calculation aie 
applied here 

• (A) Fixed weight method 

• (B) Continuously updated weights method 
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local search 


Figuie 4 4 Finding better solutions through local search method 


-V biief desciiption of both the methods is given here 


Case A Fixed weight method 

1 li(' weights aie calculated by using the fitness values of the individuals and using a linear mapping 
to hiul out the weights First, the extremes of the front on all the functions is sought and the fitness 
liuu (ion value's aie scaled using equation 4 3 Following formula is then used to calculate the weights 
ol the diffeient individuals 


(.max 


' - /»(*))/(/” 


wheie 

Wj Weight for the function 

jmax Scaled maximum value of the function m the population 
/j(x) Scaled value of the function 

Scaled minimum value of the function m the population 
M Number of different functions 

In this calculation of the weights, the objective functions are of minimization type When a solution 
IS (lose to /J”*”, the numerator becomes one, causing a large value of the weight for that function And 
(hat function is given a lot of importance Figure 4 4 shows the embedded concept If the Pareto- 
optimal front does not have the bias toward any particular region of the front, this strategy works fine 
and approximates the weights fairly well But if the Pareto-optimal front is having a bias towards any 
legion the linear approximation of the weights does not perform well 
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Case B Co 7 ifn}uously updated method 

poi some piobiein which have the bias towards the some part of Pareto-optimal front, it was ob&eived 
i]^at th(' picMous wcnglit (al( Illation vStiategy does not perfonn well, as almost all the solutions ti\ to 
iiCcU to OIK' solution and foim a cluster thus taking away the diversity This is imdesnahk' 
It IS {lu(' to tlu' itvison that the weighting function given in equation 4 4 tries to map the region with 
a lincvu luiK tion although the Pareto-optimal front may have a bias towards a paiticulai legion thus 
niMiig ('(lual impoitaiue to all the solutions But this method does not give the most optimum weight 
\{'(toi ( oi ic'spoiulmg to the difteient solutions So in order to reach near to the natural weights of th(‘ 
solution on a Paic'to- optimal front a continuous updated method is proposed 

Tins ecpiation 4 4 is based on the assumption that the solutions at the extremities are true extrem- 
ities But it may be a tiue extremity under some constraints, not a free extremity This fact is taken 
iMC while hnding the weights of individuals other than extremities, in the continuous updated weight 
(al( Illation method This method is more clear with the help of the Figure 4 5 for two objectives Foi 
inoie than two objcH fives a similar procedure can be adopted easily 



Figure 4 5 The continuously-updated weight approach is illustrated 

First, the extreme solutions A and B are assigned a weight vector according to equation 4 4 A local 
search procedure is used to find the corresponding optimum solutions Based on the new values of 
and values, an intermediate solution I is assigned a weight vector Preferably, this intermediate 
solution I can be chosen as the one having maximally away from the extreme solutions The local search 
I)io(eduie can be used with this intermediate solution and the optimum solution can be found Now 
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(vich of the two intervals ar€ considered separately region between A and I and the region between I 
and A An intermediate solution (preferably the solution in the middle of the range) is chosen in each 
th(' solution L is diosen in the range between A and I The pseudo-weight vc'ctoi foi this 
solution lakulated as follows 




XI) 






- f 


/t-/r 


(4 5) 


ulu'K' solutions a and b aie extieme solutions foi the range under consideration Thus foi solution L 
solution u IS A and solution b is I The paiameters and are pseudo-weights associated with the 
('\ti('iiie solutions Whoii pseudo-weights aie calculated for objective functions, they can be noiniali/.ed 
This pio( edmc’ can be repeated till the local search is applied to all solutions This update of pseudo- 
wf'iglits vitli a local seaich proceduie is advantageous in three ways - 


1 The tiue optimum solution corresponding to a weight vector can be found, and 

2 The calc ulated weights for each optimum solution is close to their true weights 

3 This givc's suffic lent space to spread the solutions, which are biased towards one, region i e tlic' 
PaK'to-optuiial liont is learranged according to the bias towards some region 


The dc'tailed methodology is given here 


1 Take' the* extiemities of any function and assign the weights as giving importance to that function 
only and igiioimg others Assign weight value unity to the corresponding best function value 
Use the* lact u)t = 1 

2 Find new extremities after hill climbing 

3 Set weight of the extremities as bounds on the weights 

4 To find the approximate middle point of the region bounded by the weights, weight of all the 
individual between the region is find by using formula 


w: 


W2j + 


(Wlj - W2j) * [h] - fj) 

[ft] ~ fi]) 


wheie 

'-"j weight of the it/i individual 

W 2 j Lower weight bound to the jth fitness function 

u>ij Upper weight bound to the jth fitness function 

/ 2 j Scaled value of the jth objective function corresponding to weight W 2 j 

fij Scaled value of the jth objective function corresponding to weight wij 

f ' Scaled value of jth objective function of the ith individual 
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5 Do hill climbing on t.ho individual which is approximately in the middle of the given bounds and 
find the pseudo-weight tv using equation 4 5 

(, \c)\v split tlu' H'gioii into two iiaits ' One between fiist bound and the middle point and othei 
Ix'tween the middle point <uul second bound Goto Step 4 and do the hill climbing lecuisiveh 
until all till' nidividiuils <u(' taken caie 

This foi inula of wcnght ( <ih illation ensuies that the weights are normalized oi u’l = 1 

,6 5 Clustering 

Ik' nuinbei of solutions obtanu'd just attei the iion-domination search on the solutions obtained b\ 
(al seaicli may b(' moie than lequned To reduce the number of solutions to the desiied numbei 
u' ( lusteimg appioadi is used The appioach used for the clustering is similar to that used in [36] for 
‘(hieing the size of non-dommated set The Figure 4 6 is an illustrative sketch of the clustering Here 
iitially all N solutions aio initially assumed to belong to a separate cluster Thereafter the centroid 
all the diftcieiit cliisteis and then oiuhdean distance dc between all these different clustei pairs is 
ill ulatc'd Two difleient ( lusteis having the minimum distance are merged to form the bigger clustei 
Ins jncKeduie is lejieatc'd until the numbei of different clusters is reduced to the numbei of solutions 
^cjuiicd by the us('i Finally, fioiii each cluster the solution closest to the centroid is selected Rest all 
le solutions aie diopped This is how the cardinality of the solution set is reduced Figure 4 6 shows 
K' non-dommat('d solution set m the open boxes and the reduced set m solid boxes 


f2 



Figure 4 6 The clustering method to reduce the non-dommated solutions set 


*The region is split according to the function values before hill climbing Because after hill climbing the individual 
'Ic'cted may move to one side and include or exclude few individuals whicdi were previously inside the region of one bound 
ud this approx middle point but now moved to the middle pomt and other bound This shifting from one side to other 
lay cause some error in finding natural weights 
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4.7 Closure 


\hv innlt t iv(' optimization problems have more than one objective which are to be satisfied 
sunult anoouslv d he piohlom becomes complex if the objectives are conflicting in nature The solution 
lo t li('s(' piol)l(mis IS a sot of solutions called Pareto-optimal sets Different multi-objective optimization 
algoiithnis ai(' sugge^steni a specific MOGA — elitist non dominated sorting genetic algorithm (NSGA- 
1 1 ) IS s('l('{t('d TIh' (loinmance definition used for ranking the population is changed to handle the 
( onsiiamls Ea( h solution of the Pareto-optimal set obtained by NSGA-II undeigo a hill climbing local 
s{'ai( h and the' now solution is obtained This new solution set gives a new Paieto- optimal fiont aftei 
t ht' non (lounuat ion sevue h is earned out for eliminating the dominated solutions It is expected to ha\e 
a l)('tt(u (oimngtmte than bcfioie The size of the non-dommated solution set is also expected to be 
hhIikchI as a lew nou-dommated solutions can converge to the same optima The local seaich uses a 
smgl(' ol)](‘( five function The multiple objectives are converted into a single objective by using eithei 
lh(' fixc'd weight stiategy oi by using the continuous updated weight strategy The numbei of solutions 
IS iuitlKU hhUkchI to a usei specified number by the clustering method 
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Chapter 5 


Multi-objective problems and the 
results 


In tlu' thud (haptoi, the applicability of the hybrid approach to the solve the single objective shape 
(h'sigu juobleiiis has been discussed It is found that the solutions obtained by consideiiiig single 
ol)](>(tiv(' at a tunc is not sufficient, as with a change in random number the optimum function value 
(!o('s not change much, but some times the appearance changes very much This also bung changes in 
^1 1 ('sses , 1 b well In this chapter the study is extended to the most important part of the study application 
of hvbiid appioach to the multiple objectives problems It means same problem is formulated as Multi- 
ob]V(fnH optimization problem by bringing more objectives To solve this set of pioblems, from the 
s<'t of diffeient multi-objective genetic algorithms (MOGA), Non-dominated sorting genetic algorithm 
(NSGA-II) IS selected This algorithm can maintain diversity on the front This feature is exploited to 
hud a number of diverse solutions The given solution set is then refined using a local seaich to find 
1 1 k' ojitimal solution On top of this, a clustering algorithm is used to reduce the number of available 
solutions to a desired number This helps the designer find as many different solutions as he wants and 
I hen selection can be done according to his requirement The use of genetic algorithms help the designer 
111 finding new shape which increases the creativity in the designs This has the biggest advantage of 
giving the multiple results in just one shot 

For this study all the problems which were considered for single objective case are solved for two 
objectives The problem of design of cross-section does not involve any loading so the weight of the 
design does not affect the design All these problems other than the problem of design of cross-section, 
aie solved for two different cases 

1 When the weight of the design is not considered 

2 When the weight of the design is also taken into consideration while calculation of the loading 

For the problem of the design of the cross-section, moment of inertia to weight ratio for two axis, is 
maximized by assigning one objective function to each axis All other problems solutions are solved for 
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lollowing two objoctivph 


• Fust objoctivp IS nnnimization of weight 

. S<>< und ub,c.c t.ve is mumm/atiou of 1/Stiffness Since stiffness is given as inveiso of displacement 
St) tins objet tiv(> ( <ui be ledohned as minimization of displacements 

Foi all the- test piobk-nis c onsideied, other than the design of cross-section, following GA pararneteis 
ait- selet tt-d 


l’o|nil<it ion Si/t- 
C'iosso\(‘i Piobability 
Miit.itioii Piobability 
\o of Geneiations 


30 
0 95 


1 /String length 
150 


loi the- design of cioss-seclion the only difference is that the population size is taken to be 100 instead 
of the 30 

The- dc'sciiption of the problems is given in the third chapter The change in the solution is that 
t,h(> maxiinum displacement is also treated as an objective here The deflection is scaled such that the 
oidc'i of both objective values becomes the same This objective is also a constrained objective 

5.1 Design of the cross-section 

The pioblem is described m the chapter 3 The results for the single objective optimization are also 
jnesented in the chapter 3 The results obtained for these problems for the single objective case clearly 
show that the approach considered has the potential of finding the optimal solution The next and more 
important step in solving the problem is to consider more than one objectives This problem is solved 
foi two objectives- maximization of moment of inertia to the weight ratio for two axis 

First plot in the Figure 5 1 shows all the non-dominated solutions obtained after the NSGA-II 
simulation runs Since the population size is 100, almost 100 non-dommated solutions are obtained 
Second plot shows the solutions obtained after the local search is earned out on the solutions set obtained 
aftei the NSGA-II simulation runs This is clear from the plot that the cardinality of the solution is 
1 educed The non-dommation local search on this set of the solutions gives the final non-dommated 
set as presented in the third plot The number of the solutions present in this set is 81 The graph 
shows that the diversity is lost on the Pareto-optimal front The solutions which are maximally away 
fiom the ends, mostly give way to the solutions near to the optima for one objective functions Here 
some solutions have converged to one particular solution after the local sehrch because the weight for 
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(vuh solution IS detci mined individually, without considering the effect of the others That is also the 
KMsoii why soiiK dominated solutions are obtained after the local search Out of these 81 solutions 
iiinc solut ions ai(> sel(>< tod by means of the clustering These are shown in the last plot That’s hov the 

iiihU si'ltH lion of tlu' solutions is made 



Scaled MI U) weight ratio lor x-axis Scaled MI to weight ratio 1 or x-axi*! 


Figmo 5 1 Hybrid procedure to find nine trade-off solutions for the cross-section design problem 

In oi dei to visualize the obtained set of nine solutions having a wide range of trade-offs in the scaled 
moment of inertia (MI) to weight ratio for both the axis, the shapes m Figure 5 2 are presented It is 
< leal that starting from a solution with large MI to weight ratio for one axis, how the shapes for the 
laig(' MI to weight ratio for other axis are obtained by the use of the hybrid approach The solutions 
<u(' jnosented in the order of increasing MI to weight ratio for one axis, m the form of a 3 x 3 matrix 
Tlie fn&t solution at position (1,1) in the matrix, has the maximum MI to weight ratio for x-axis The 
solution places the vertical bar on ends to increase the moment of inertia for other axis, i e increase in 
MI to weight ratio for other axis The solutions (1,2) and (1,3) increases the ratio for y-axis by putting 
more material far from the center Sixth solution (2,3) of the matrix, gives more importance to the 
moment of inertia to weight ratio for the y-axis The solution has lesser moment of inertia to weight 
latio for x-axis The ninth solution (3,3) is having the maximum MI to weight ratio for y-axis It is 
dear from the solutions that it tries to keep more material far from the axis to increase the desired 
1 atio This also tells the method of placing material to increase the moment of mertia in a optimal way 
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□ c □ 

n n n 

Figure 5 2 Nine trade-off shapes for the design of the cross-section 
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5.2 Design of the cantilever plate 

Th(' l().Kling (oiKhtions and the supports for the cantilever plate are shown m Figure 3 1 The results 
fdi till (<is(s wluii the weight of the plate is considered to have some effect on the loading and when 
ilu- vmght oi tho d<>signod plate is not considered are presented separately 

5 2 1 Design of the plate when its weight is not considered 

I i..,ui( .) 3 pi(S( iits the fom stejis of the pioposed hybrid method applied on this pioblein Fust plot 
slums f li(> ii()ii-d()iiimat('cl solutions obtained using NSGA-II Since the population size is 30 NSGA-II 
IS able to find 30 dificicuit uoii-dominated solutions Thereafter, the local search method is applied 
lioii! (vuh iion-doiiimatod solution and new and improved set of solutions aie obtained The thud 
plot IS th(' K'siilt of th(' iioii-doiiimated check of the solutions obtained after the local search method 
1 111 ('(' doiiiiiiated solutions ait' eliminated by this process The final plot is obtained after the clustering 
opi'iafion with a (hone of nine solutions The plot shows how nine well distributed set of solutions 
ai(' found hoiii the thud plot of 27 solutions If fewer than nine solutions are desired, the clusteiing 
iiK'chanism <au Ix' set accoidmgly 

111 oid('i fo viMialr/e the obtained set of nine solutions having a wide range of trade-offs m the 
weight and scaled defic'ctioii values, the shapes in Figure 5 4 are shown It is clear that starting from 
a low-w('ighf solution (with large deflection), how large-weight (with small deflection) shapes are found 
by the hybiid method 

It IS inteiesting to note that the minimum weight solution eliminated one complete row (the bottom- 
most, low) in Older to reduce the overall weight The second solution (the element (1,2) in the above 3x3 
mat 1 ix) ( oiK'sponds to the second-best weight solution It is well known that for an end load cantilever 
plate, a parabolic shape is optimal Both shapes (elements (1,1) and (1,2)) exhibits a similar shape 
As the importance of deflection increases, the shapes tend to have more and more elements, thereby 
making the plate rigid enough to have smaller deflection In the middle, the development of vertical- 
stiffener is interesting This is a compromise between the 110111111111111 weight solution and a minimum 
deflection solution By adding a stiffener the weight of the structure does not increase much, whereas 
the stiffness of plate inci eases (hence the deflection reduces) Finally, the complete plate with right top 
and bottom ends chopped ofif is the mimmum deflection solution 

All nine solutions (and if needed, more can also be obtained) with interesting trade-offs between 

1 

weight and deflection are obtained in one simulation run of the hybrid method More importantly the 
solutions to this problem present the way of putting material to reduce the deflection in an optimal 
way 





5 2 2 Design of cantilever plate when weight is also considered 

\()\\ t ht' of the designed plate is also considered to effect the loading Figure 5 5 shows the 

iion-donnihitt'd solutions obtained using the hybrid appioach A caieful look on the solutions foi the 
( as(' wIkmi wc'ight is (onsideied and when weight is not considered clearly tells the effect of the weight 
{ oiisulc'i at ion I lu" high density increases the load so the solution spread is reduced on both axis i e the 
iaug(' of tlu' w(aght and the lange of the deflection is reduced The NSGA-II gives 30 non-dominated 
soliUions riu' local seaich takes the solutions to a better distribution The local seaich impioves 
ili(' di\(‘isit\ among the solutions The number of the solutions is leduced to 20 fiom 30 aftci the 
iioii-doiuniatiou (heck is carried out on the set of the solutions obtained after local search Clustenng 
nH'(h<misni is used to extiact the nine diversified solutions from this non-dominated set These nine 
solutions ai(' shown in the Figuie 5 6 

! h(' solut ions aie aiianged in the order of the increasing weight The minimum weight solution (left 
solution 111 the top most low), exhibits the shape veiy similar to the parabolic shape The solution 
has two (hcnnents ]ommg the node bearing the load to the lest of the plate Next solution element 
{ I 2 ) has iiioio symmetric shape Another interesting shape is the third solution (showp m the (1 3)th 
position of the' matiix) This solution has a vertical stiffener besides the parabolic shape This gives it 
<ulditu)nfd stiemgth without increasing the weight much Other solutions increases the mateiial iieai to 
lh(' stiflc‘U(‘i to me lease the strength and hence reducing the deflection Finally the solution piesented 
in the' position (3,3) of the matrix is for the maximum weight Here the number of stiffeners has also 
me Kvisc'cl so the deflection is reduced This solution is significantly different than the maximum weight 
solution ohtamed m the previous case Most interesting point of this solution is the number of holes 
riic' numbci of holes is increased due to more number of stiffeners This help in reducing the weights 
while keeping the deflection minimum 




Figure 5 4 Nine trade-off shapes fee the cantilever plate design 
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NSGA-1 1 
Non-dominated local search 
Clustered Solution 


20 25 30 35 40 45 50 

Weight 


Fignio 5 5 Hybiid procedure to find nine trade-off solutions for the cantilever plate design piobleni 
w Ih'u weight IS also considered 


Figure 5 6 Nine trade-off shapes for the cantilever plate design when weight is considered 
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5.3 Design of the simple supported plate when point load is applied 
on the topside 

\('\t piohlc in IS the ck'sign ot the simple supported plate when a point load is applied on the topside 
lit the pl.iti as slumu lu the Figuio 3 2 The results for the cases when the load of the design itself is 
.uiouuied ioi the (aleulation of the weight, and the case when load is not accounted as loading, aie 

jiK'st'Uted 

5 3 1 When weight of the designed plate is not considered 

[■ignie a 7 shows the hybiid appioach to find the nine solutions The NSGA-II obtains the non- 
doiiunati'd sot of the solutions These solutions are further improved by the local search giving a 
hot (c'l distiilmtioii on the Pareto-optimal front The shape of the Pareto-optimal front is also improved 
HI t his (>x('K ISO of cairyuig out the local search The number of solutions is reduced to 22 after carrying 
out a nou-doiinnation check on the solution set obtained by the local search The method reduces the 
( ai(hiiaht.v of the pioblem to some extent Nine solutions are selected from this set using the clusteiing 
appioach These solutions aie presented in the Figure 5 8 

An analysis of the solution leveals some interesting information The minimum weight solution 
( 1 1 )th 111 the matiix, tends to use one row less but since the load is applied in the top row, so in-ordei 
lo caiiy the load one element is placed m the top row The third solution is an interesting solution 
A < <ueful look of Figure 5 7 reveals that this solution is a knee solution It means the solution makes 
,i huge sacrifice in the deflection to achieve a small advantage m the weight Similarly, to achieve a 
Miiall advantage m deflection-loss, a large sacrifice m weight is needed Shapes m position (1,2) and 
(2.1) can be compared with respect to the shape in position (1,3) Shape m position (3,1) or solution 7 
IS also interesting In order to have further reduction in deflection stiffening of the two slanted arms is 
needed These stiffener of the slanted arms form a close loop or it can be said that the solution has got 
i holes Solution (3,2) or 8th solution has thickened the slanted aims, but one hole is there to reduce 
tlu' weight Finally, the absolute minimum deflection shape is the complete rectangle with maximum 
possible weight One very important point to note here is the almost symmetric nature of the shapes 
Starting with the minimum weight design having two slim slanted legs down to thickening the 
legs to make them stiff, followed by joining the legs with a stiffener, and finally finding the complete 
1 octangular plate having minimum deflection are all intuitive trade-off solutions In the absence of any 
sue h knowledge, it is interesting how the hybrid procedure with NSGA-II is able to find the whole family 


of different trade-off solutions 



5 3 2 Weight of the plate is also considered. 

Njow di<' plfiU* IS designed when the weight of the design is also considered while calculation of 

tlu' stu'ssf's and deiiei tioii Pigme 5 9 shows the solutions obtained by the NSGA-II simulations Local 
si'.uili on lliesc solutions iiu leases the diversity The non-dommation check gives 29 non-dominated 
siilulums iillei the !o( al seaich Here the number of solutions is not leduced veiy much but the spread 
()1 the solutions is impioved The solution set is reduced to nine by clustering The final nine solutions 
.lie pieseiitc'd in the tonu of a matrix of 3 x 3 m Figure 5 10 

A ( lose' look on tlu' Figuie 5 7 and Figure 5 9 reveals the fact that the diversity of the solution set 
IS lediiced whc'ii tli<' wc'ight is also taken into consideration Now the spread of the front is leduced 
to 13- H) units mstcvid ot 15-60 units on the weight axis and the spiead is also reduced on the scaled 
delic'c t ion axis I’lic' shape of the Pareto-optimal front is evident of the improvement in the solutions 
hv the’ Inc.vl sc'aich The solution in the first rov/ is the same as the solutions when the weight is not 
(oiisidoied The solutions in the second row le solutions (2,1) and (2,2) put more material on the 
aims Since' the' weight of the design is also taken as the distributed loading so the material is also 
put in tlu' optmwl way such that the effect is not to increase the deflection The material is put near 
to the' supports siuli that it has lesser effect due to the weight of the designed plate The solution at 
the' iiositioii (2,3) soc'm very interesting as this solution puts same amount of the material everywhere 
Re'st, all other solutions in the last row of the shape matrix seem a derivative of this solution All the 
solutions 111 the last row have lesser deflection and more weight, but unlike the case when weight is not 
( onsidered, the designs put lesser material m the design and the reason is very obvious Here more 
uiatc'iial means moie stress and more deflection, hence lesser material is used So the solution with 
maximum amount of the material is not the same as the whole plate The solution with the maximum 
weight IS also now having the deep chopped comers with thick arms Interestingly, all the solutions are 
symmetric in nature and no hole is present m any solution 


68 


NSGA-n-^ 
Non-dominated local seaich D 
Clustered Solution ♦ 


15 20 25 30 35 40 45 50 55 60 
Weight 


ignu' 5 7 Hybnd procedure to find nine trade-off solutions for the simple supported plate design 
ohk'iii will'll wi'ight IS not considered 


Figure 5 8 Nme trade-off shapes for the design of simple supported plate with a pomt load apphed on 







Figure 5 10 Nine trade-off shapes for the design of simple supported plate with a point load applied 
cm the top and the weight is design also considered 
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5.4 Design of the simple supported plate with a distributed load 
applied on the top 


t 


i ho snpiHiits and t!i<' loatlmp, conditions of the problem are shown m Figure 3 3 The pioblem is such 
hai If loquiios ,ill tlu ( k me nts in the first row to be present to make the design geometiically feasible^ 


1 ho u'sulfs ii'i liofh of the cases aie presented in the next subsections 


5 4 1 Weight of the plate is not considered 

hast tlio losnlts foi tlKMasc", when weight of the design is not consideied, are piesented m Figiiie 5 11 
rh(' solutions obt.unod by the NSGA-II are presented here The Pareto-optimal fiont obtained Iw 
\SGA-II luiis li.nc' a step ui the front Local search improves the spread of the solution set and 
iho (oii\o\ Pai('to-opf luial front is obtained after the local seaich The spread on the weight front 
IS iiiipioM'd fioin 26-53 units to 22-60 units and on the scaled deflection the impiovement is clcaih 
Msible as tlu' lesult of local search The non-dommation check on this new solution set, yield 21 well 
spac ed iion-doimuated solutions This reduces the cardinality of the problem and take the appioach 
t()W<u(ls uioK’ juartical side Nine diverse solution are selected through clustering and are presented iii 
flic' FiguH' 5 12 

'Phc' solution with the least weight is solution (1,1) of the matrix This solution has the straight 
ana )omiug tlu' l(*ft end support with the top load carrying part of the plate The right arm has got a 
small twist at the joining point The second solution at (l,2)th position is a symmetric solution which 
shifts th(‘ joiiiiug points of the supports with the load carrying arm towards inside giving a shape of the 
l<ible From the solutions shown m Figure 5 11, it is found that this solution is a knee solution A big 
(ompioimsc is made on the deflection, to accommodate more weight The solution (1,3) and (2 1) keep 
t hc' symmetiy but thicken the arms The trend of keepmg more material on the supports continues such 
that the tabular structure is obtained and the supports are joined with the top row elements through 
the straight arms The small shift m the symmetry towards one side is just due to the attempt of the 
scvuch of the solutions with different stiffness The maximum weight solution is the whole plate (the 
iiiaxiimim possible weight) 


5.4.2 Weight of the designed plate is considered 

Figure 5 13 shows the result obtained by using the hybrid approach to solve the problem Initially. 
NSGA-II is used to find the solution set of the problem The nature of the Pareto-optimal Ftont 
obtained by NSGA-II gives an expression that the some part of front is concave The solutions are 
improved by the local search The Pareto-optimal front is improved by the local search and the front 
i.s found to be convex with the presence of a few steps Only 17 non-dommated solutions are obtained 
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1)\ this pioce.sh of local search and the non-domination search Actual number of different solutions 
IS {ouiid to be nine Thus the number of solutions is reduced fairly The solution set is less diverse 
ih<ui the' ( as(' when weight of the design was not taken care The weight axis diversity has leduced to 
iO-oO lioin 20-bO The diveisity on the deflection axis has undergone a major change as it is reduced to 
1 i)-,l 2 fioin 1 0-8 0 in the pievious case The reduction in the number of different solutions may be the 
pu'sc'iH (' of optimal points in the basin wheie the different solutions are taken by the local search Thus 
no ( oiitiol on spread ( an be obtained by this weight selection strategy Besides this the wc'ights for the' 
-,olutious m the (oiitave region may not be very appropriate This is the limitation of the weighted sum 
stiati'gv of not reaching the optimal weights m the concave region The solutions which are very near 
lo ('ach other also converge to one solution That’s why the diversity after the local search is affected 
\me diffeient solutions are presented m Figure 5 14 

Tlu' mimmuin weight solution (position (1,1) solution) is symmetric and has the got the tabular 
sh.ipi' Other solutions like those shown m (1,2) and (1,3) position of the 3x3 matrix have more 
matc'iial placc'd on the aims The thud solution is a knee solution as a small change in the weight has 
( ausc'd a large change in the deflection values The solutions in the second row are veij near to each 
otheu 111 the fitness space and therefore have very much similarity m the shapes Solutions in the last 
low restore diveisity in the fitness values and as well seen in the shapes More weight is put on the 
ai ms sue h that the joints are thickened The tendency is to keep the weight or the material near to the 
Mijipoils thus reducing the effect of the increase m weight This way the solutions does not have got 
\('iy much increase in the deflection The maximum weight solution have thick symmetrically placed 
M ms All the solutions are almost symmetric and have a near tabular shape with thick slanted arms 
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Weight 


Fi^uu' f) 11 Hvbn<l pioa^duie to find nine trade-off solutions for the simple supported plate design 
])ii)bkuii witii <i distnbuted load on the top and the weight of the design is not considered 



Figme 5 12 Nme tradeoff shapes for the design of a snnple supported plate witha distnbuted load 

applied on the top 
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FiEuie 5 14 Nme tiadMff shapes for the desigi. of sunple supported plate with a distributed load on 
the bottom side when the weight of the design is also considered. 


74 




5.5 Design of the simple supported plate when a point load is applied 
on the bottom side 

riu' loaduiR <uh 1 till' suppoits toi this problem is shown in the Figuie 3 4 Now the loading is applied 
on the lowc'i si(l<' of the plate, which makes it feasible to eliminate the top portion, if not required The 
s.utu' tK'iid is ()hs('i\('(l ill the solutions for both of the cases, when weight of the design is considered 
,ui<l when weight of IIk' design is not considered 

j 5 1 Weight of the designed plate is not considered 

1 h(' lesulf s toi t Ins jnohleiii is shown in Figure 5 15 Initial solution set of 30 non-dominated solutions 
IS olitaiiK'd b\ tiu' NSGA-II Foi low weight the solution set obtained is having the front which seem to 
hav(' soiiK' st('])s Tlu' solution set has a few cluster of the solutions The local search on this solution 
si't tuithei inipiovc's the diversity and the shape of the Pareto-optimal front is also found to be convex 
w it bout steps D('spit(' the clusters of the solutions in the initial solution set, the new set after the local 
si'anli has sulutaon s]neatl at intervals on all over the front and 28 non dominated solutions are obtained 
aftei the' non doininataon check Out of these 28 solutions only 13 are different solutions The spread 
of t,lK' solut ions on t h(' weight axis is from 10 units to 60 units Similarly, on the deflection axis, the 
spu'.id IS fiom 2 0 to if 0 units Nine different solutions axe taken from the set of 13 different solutions 
and piesenU'd in Figme 5 16 

Tlu' solutions obtained here are very interesting The minimum weight solution at position (1,1) 
111 the rnatiix eliminates the thiee rows completely The solution is very interesting in the sense that 
I his joins the (uids in such a manner that there are two holes These holes in the structure reduce the 
wc'ight Tli(' thud solution is also interesting The solution is symmetric and is placed such that the 
matoiial is jilaced inoie m the middle of the plate where the load is applied This solution follows the 
shape of the bending moment diagram for this problem More material is placed in the zone of higher 
bending moment The solution is also interesting in the manner that this solution reqmres very much 
1 eduction m the deflection to optimally increase the weight by a small amount as shown by the slope 
of the line joining this solution and the previous solution The fourth solution at position (2,1) has 
jnit more material m the middle of the plate near to the load This helps in reducing the stress by the 
increase in the cross-section area on this section Similarly, the fifth solution (2,2) is also symmetric 
and have moie material near the load The arms are also thickened now This solution has made use 
of one more row but the material is placed symmetrically Interestingly, all these solutions maintain 
two cavities and try to maintain the symmetricity in the shapes The solution for the maximum weight 
have whole of the plate as the solution 
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5 5 2 Weight of the designed plate is considered for the calculation of the loading 

\()\\ t ho < as(' when the weight of the designed plate is also considered for the calculation of the stresses 
uid (iolh'c tioii is <inalysed The lesults are presented in the Figure 5 17 The solution set obtained b\ 

1 Ik' \S(i V-11 gives loasonable diveisity but the shape of Pareto-optimal front is not smooth The local 
st<n(li unpiovt's the solutions such that the resulting solution set gives a Pareto-optimal set without 
^K'ps aiui tlic' (liveisity on the front is impioved 29 non dominated solutions aie obtained fiom the 
IK ni-doinnhil loii seauh on the set obtained after the local search, out of which only 15 aie difteient m 
slia|)(' ! 1 h' spuvul of the solution set is slightly lesser than the case when the weight of the designed 
piat(' is not ( onsKhnod foi calculating the load, as the loading increases with the increase in the weight 
solution which are sufficiently distinct in the objective function space are selected and piesented 
111 t h(‘ Figinc' 5 18 

riic' unniniuni weight solution is almost the same as m the case, when the weight of the plate is not 
a< < ountc'd foi c alculation of the loading Solution (1,2) and (1,3) put more material in the middle of the 
platc' But thc'y aie no more symmetric The solutions has a tendency to increase the cioss-sectiou ]usi 
.il)o\ t !i(‘ load to reduce the stresses and the deflection The fourth solution (2,1) seem yeiy mteiestmg 
.IS this lias put the maximum material in the middle just above the load The material experience a 
pull tuw.uds down m the middle and the slanted arms restricts the extent of the pull Othei solutions 
ai{‘ sue h that they use more mateiial to stiffen the slanted arms This make it possible to leduce the 
dc'fic'c t ion But the cavities on the bottom side are retained This shows that the material at the place of 
( avit K's ( an be safely taken away to reduce the weight The maximum weight solution is also having the 
( <ivd K's and remove the material on the top ends Interestingly, all the solutions are nearly symmetric 
in natuH' and the material is placed according to the variation m the bending moment 
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Weight 


l’i<j,uH' 5 15 Hybud procedure to find nine trade-off solutions for the simple supported plate design 
piohk'iii when point load is applied on the bottom side and weight of the design is not considered 



Figure 5 16 Nine shapes for the des«n of snnple supported plate with a pomt load applied 

on the bottom side 
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Figure 5 18 Nine trswic-off shapes fox t]i6 design of simple supported plate with a point load applied 
oil the bottom side when weigt of the design is also considered 
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5.6 Design of the simple supported plate with a distributed load 
applied on the bottom 


1 1„. I, .m.l snpp,n . conditions tor tins problem is shown m Pignre 3 5 The loading is applied such 

iliut 111.' ..II cl....... Ills ... Iho last low must be present to make the solution feasible The solution lo. 

, lus piobleii, 1 .... elm....,itt. the lows from the top rf needed, without making the solution geometricallv 

inicvisihle 


5 6 1 Weight of the plate is not considered 

I'nsi llu' pioblfMn IS solvent! foi the case when the weight of the designed plate is not taken, while 
( al( ulatnig th(> stiessi's The lesult for such case is presented in Figure 5 19 The solution set obtained 
fioiii th(' NhGA-II 1 mis is laiily diverse and the shape of the Pareto-optimal front is also nearly smooth 
riu' defiec tioii toi this pioblem is found to be more than other problems These solutions undergo hill 
( limbing local seaidi and the solution set is further improved The solutions after the local search have 
gi)t th(' cuds of the Pan'to-optiinal front stretched by a small amount The rion-dominated set of 25 
solutions foiiuh a p('ii('( t Pareto-optimal front Effective number of the different solutions is lesser than 
I h(' iioii-dommati'd solutions as many solutions from the NSGA-II runs have converged to one solution 
C'lustcuiig helps to eliminate the solutions which are very near m the objective space or are co-incident 
S(‘t ol nine difitumit solutions found by clustering is presented in the Figure 5 20 

The solution set obtained is very interesting The solution with the minimum weight have the last 
low only R('st whole ol the plate is eliminated This solution is the same as found for the single 
obji'ctive cds{' solution while optimizing for the minimum weight The second solution, element (1,2) ol 
th(' 3x3 matiix, has more weight placed symmetrically This follows the shape of the bending moment 
diagram to have an almost uniform distribution of the stresses This solution places the material in a 
( (impact manner such that the deflection is reduced substantially Next solution (1,3) also exhibits the 
same trend But the improvement m the weight and deflection is not as much as it has been in the 
s('f Olid solution Solution (2,1) and (2,3) are having a cavity in the middle This helps in reducing the 
weight without affecting the deflection much The intermediate solution is also present m fifth solution 
(2,2) The solution (3,1), (3,2) and (3,3) exhibit much difference in the weights but the difference in 
deflections is not much This shows that putting very much weight in the solution will not help much 
to 1 educe the deflection beyond a hmit The solution with the maximum weight is having the complete 
plate as the solution In this design, almost all the solutions are symmetric and the spread obtained is 
also very good 
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5 6 2 Weight of the plate is also considered 


\o\\ th<' w(nght of the design also increases the loading as the weight of the design increases The lesuh 
1(11 1 Ins ( ,is(‘ IS pit'sented in the Figure 5 21 The solution set of 30 non dominated individuals obtauiod 
1)\ ilic \SCi V-ll simulation urns is worked on by the local seaich method The local seaich helps the 
solutiiiii to gi't <1 iiioie stietchod Pareto- optimal front But on the front many solutions aic' found to 
( ( m\ ei ge to ou(' point 1 1ns leduces the number of different solutions to 15 This is due to the fa< t that 
ilu' wi'ights aie (leteiiuiiied without taking diversity of solutions into account These noii-doiiiinated 
Milutuius loim a smooth Pareto-optimal front The spread of the solutions is slightly different than the 
( ,isi when th<’ wi'ight is not considered This is reduced on the weight axis and more interestmglv the 
1 iiige on the defiiH tioii <ixis is increased Nine of these solutions aie taken and aie piesented in th<- 

1 igine ■) 22 

1 Im' solution with the minimum weight is similar to the case, when weight of the plate is not 
( ( insuieu'd This eliminates four rows from the top Infact, all the solutions in the first row of the 3x3 
inatiix <ue the same as the previous case when weight is not accounted for Solution (1,2) and (1 3) 
1( illow f h(' saiiK' tioiid and have the symmetric solutions The second solution is such that there is iniu h 
H'diu lion in the deflection as compared to the weight as evident from the Figure 5 22 Solution (2 2) 
aii<l solution (2 3) have got one hole This hole help m reducing the weight Though these solutions 
,u(' not svinmetiic the hole causes the change in stresses more, as this is experiencing a umdiiectional 
load Mok' luteiostmgly, the solutions (3,1) and (3,2) have more number of cavities It will be a matter 
ol luteu'st t ) find out the effect of the multiple holes on the solution apart from the reduction iii the 
weight of the solution The solution with the maximum weight is having a cavity and also it has got 
some mateiial cut on the corners, making the solution different from the maximum weight solution for 
them <vse when the weight is not considered 
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>MU(' 5 19 pioceduic to find nine trade-off solutions for the simple supported plate design 

oblcHU whc'u the distiibiited load is applied on the bottom side and weight of the design is not 

nisidc'i ('d 


j ff for the Simple supported plate design when a distributed load 

Figure 5 20 Nine trade-off shapes for the simpie &upp 

applied on the bottom side 







lMu,nH' f) 21 Hybiid proceduie to find nine trade-off solutions for the simple supported plate design 
piolili'iii with a distributed load applied on the bottom side and weight is considered 



P.,ure 5 22 Nme trade-ofl shape, for the smrple supported plate destgn when a drstrtbuted load 
applied on the bottom side 
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5.7 Design of the hoister plate. 

1 lu> lo.idmi- .iiui 11 k' suppoit conditions for this problem are shown in the Figure 3 6 The size of the- 
t.ikcn IS (SI) - ()() min This is a tough problem so some problem information regarding the cavity 
u.,s n<M essai N ami hem c specified in the problem This is the problem, where the strategy of having 
lixed weights has not hemi found working, hence the strategy of continuously updated weight is applied 
li('i(' and Iniiiid to woik well The solution for both the cases when the weight of the design is consider c'd 
and when tlu' weight is not considered are presented heie 

5 7 1 When weight is not considered 

kigiiii' a 21 pu'smits set of 30 non dominated solutions obtained by the NSGA-II runs The local 
siviu h 1)11 fins s('i ol solutions gives 21 non dominated solutions The solution set after the local search 
IS K'diued to a small iiumbei and further more the number of different solutions is lesser than this 
minih('i rh(‘ solutions are such that there is a sharp bend m the Pareto-optimal front This solution 
set isiedmc'd to nine solutions which are different in objective space as well as are diverse The different 
tiad('-otl solutions are shown m Figure 5 24 

riiis IS a vmy very tough problem as we are trying to get a shape with a structure open at one 
('lid whu h gi\(' iis(' to unbalanced moments, thus causing higher stresses and deflection The solution 
s('t obt aim'd fiiuilly shows the ability of the method to solve this problem also The solution with 
the niminimii weight (solution (1,1) of the matrix) obtains the shape of a hook, very often used foi 
llu' jmipose of tlu' hoister Most interesting thing about this solution is the opening on the left hand 
side without spc'cifying any problem information of the same It is also interesting to note here that 
tlu' sohitaoii foi tlic' minimum weight does not have the leftmost and rightmost columns and bottom 
most low Second solution is also interesting as this reduces the weight by a small amount but the 
change in the deflection is very much This is due to the change in the configuration The new shape 
IS hkc' the gi ippei of some arm This is connected to the upper side of the plate, hence there is not so 
time h unbalanced moment as is experienced m the first solution Third solution reduces any unbalance 
due to the loading by creating a symmetry of loading This solution further reduces the deflection 
witdiout affecting the weight by a large amount Rest all other solutions do not experience big changes 
111 deflection despite of large changes m the weight of the design Solution (3,1) is a symmetric solution 
This solution smartly eliminates the part of the left and right columns and also the corners on the 
bottom aie chopped The solution with maximum weight takes the maximum amount of plate that 
has Ix'oii allowed The last three solutions have got a big change in the weights but the change in the 
dc'flc'claon is not too much, showing the places where the weight is present without contributing much 
t,() the problem This information helps m optimizing the solutions for unknown cases 
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5 7 2 Design of the bolster when its own weight is also considered 

I ho solnl U)11 sot of 30 non dominated solution obtained by NSGAdI runs is presented m Figure 5 25 The 
p,u<'lo-u|>t niial lu)nl has a step Interestingly, only 22 of these non dominated solutions aie difteient 
Hut llu' <ii\tnsit\ on the fioiit is pietty well The local search gives 23 non dominated solutions out ot 
whuh i ht'K' OH' onh 11 distinct solutions Nine solutions are taken out of the local set and aie piesentcd 
III I li(' l''i<;ui<' .■ 2 b 

1 111 soiution s('t ha.s got mteiesting shapes The solution set with the minimum weight is having 
,111 opi'iiuig nil t lu' lett hand side This has the arm joining the load carrying part to the support This 
makes ilu' solution ol minimum weight The weight of this solution is even lesser than the miniinuiu 
ui’ighi solution of the previous case The solution (1,2) is also inteiestmg as there is some reduction 
111 I he di'ilei lion with the weight increase Second and third solution make large compromise on the 
<1('11('( liou axis to accommodate small increase m the weight The solution (2,1) gets a twist in the 
uiii )oinnig the loaded part and the support hence the deflection is reduced more causing a step m 
ili( Paic'to-optimal fiont Other solutions do not have much difference in the deflection values with 
the wc'iglit change Inteiestingly, many of the solutions has eliminated the last row of the consider ed 
scan h sp.uc' It is mteiesting to note that the solutions which puts material in the last row docs not 
liave iiiuc h dific'K'iic (' ill the deflection values Thus last row material can be safely eliminated to leduc e 
the weight, without effecting the deflection significantly (if required) The solution with the maximum 
wc'ight IS again interesting as this solution does not make use of whole plate, instead it eliminates the 
solut ions on the left and right columns The solution is interesting as the increase in the weight is very 
nmc h but the coi responding change in the deflection is not much for last five solutions Thus very steep 
slope' wit li ic'spec t to one axis is obtained This solution is different than the previous case solution of 
iiiriximum solution due to the reason that the weight is affecting the deflection also in this problem 
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nil' T) 23 Hybiitl piocedure to find nine trade-off solutions for the hoister design ' pi oblem when 
ghi of th(' jilate is not considered 





Figure 5 24 Nine trade-off sha^ies for the hoister plate design 
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5 25 Hybrid procedure to find nine trade-off solutions for the hoister design pioblem when 
wc'iglil IS also considered 







5.8 Design of the bicycle frame. 

Pif^uK' 3 7 shows the loading and the supports m the frame under consideration This problem is veiv 
nitc'K'stiug as diiect application of this work is presented in the real world problems The design of 
(i aiuc IS made with a vciy high static load than usually experienced The hybrid approach of the design 
givc's a numbei of solution which gives not only optimal shapes but also the innovative ideas for the 
diftcH'iit shapes of the bicycle 

5 8 1 When weight of the bicycle frame is not considered 

iMguu' 5 27 present the results of the NSGA-II simulation runs The diversity of the solutions on the 
Pau'to-o])tnnal front is very good The non dominated solutions obtained after the local seaich form a 
new lefined Paieto-optimal front The obtained Pareto-optimal front clearly resembles the initial idea 
d('S( ribed The spread of the solutions on the new obtained Pareto- optimal front is good, and the range 
of fitness values is improved Out of 25 non-dominated solutions there are 20 different solutions Five 
sohilioiis from this widely spread front are selected and presented m the Figure 5 28 Lesser solutions 
au' (onsidcred here, to avoid the cluster of solutions available for the designer, which ngake the task of 
decision making difficult 

The shape of the bicycle frames obtained by the hybrid approach are presented m the form of 
the bicycles The shapes obtained are very interesting The solution for the minimum weight is very 
inteiesting as this has a vertical bar joining the seat with the paddle bracket and the rear wheel such as 
maximum load is taken by this part of the design The seat is joined to the handle by another bar This 
(h'sign IS not available in the market, but it can be a good solution to consider Next shape obtained is 
a del ivative of the minimum weight solution This solution thickens the vertical bar joining the seat 
and the rear wheel and the paddle bradket Shght modification in the weight reduces the deflection by 
a laige amount It is interesting to see how the hybrid approach eliminates the material to find optimal 
solutions Another solution has got a slanted bar joining the handle and the paddle support This 
solul.ion IS sinnlai to the most common cycle available in the market The thickening of the arm joining 
the tear wheel and the seat make it stiffer by reducing the stresses and hence strains Most interesting 
thing about this solution is the cavity obtained by the solution The solution with the maximum weight 
eliminate all the cavities The steps are found near to the front wheel This is due to modelling of the 
fiont wheel space as a staircase This solution has eliminated all the material in the last row but due 
to the lequiiement of paddle wheel bracket, one element is present m the last row 
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5 8 2 Design of the bicycle frame when its own weight is also considered 

I his (k'sigii IS moie near to the real life design The Figure 5 29 presents the solutions obtained by 
tlu' \SGA-II simulation runs This gives 25 different shapes of the bicycle frames The solution set 
uucl{‘ig()('s the local seaich and 15 non-dommated solutions are obtained Local search mci eases tlu' 

1 ()1 t he solutions obtained on the weight axis but the lange is considerably reduced on the deflection 

<i\is as IS evident from the Figure 5 29 Most interestingly the number of distinct solutions is 8 Tlu' 
small numbei of solutions obtained show the efficiency of this approach for reducing the cardinality ol 
this piohlem Out of these solutions five distinct solutions are selected using the clustering appioach 
( 'onsulca ation of the weight has got a significant increase m the deflection as compared with the case' 
w luMi th(' weight of the frame is not accounted These solutions are presented in the Figuic 5 30 

Ih'K' the shapes obtained are very interesting as it gives multiple solutions m one simulation lun 
UK hiding some familiar shapes The solution obtained for the minimum weight has got a vertical aim 
joining t he seat with the paddle support This paddle block is joined to the handle support and the leai 
w hcc'l biacket The similar solution is seen m the day to day life The next solution is one step ahead 
of the minimum weight solution in the sense that this thickens the joint of the verticjal arm joining 
paddlc' block and the paddle wheel with the rear wheel bracket Thus the solution becomes moie stiff 
1 h(' solutions shown in the second row are also interesting This has a huge cavity in the solution 
w it bout giving any problem information to the solution The solution is similar to one of the models 
availabk' m the market, but it eliminates the bar joining the paddle block and the seat This shows 
1 h<it tins aim can be eliminated to reduce the v/eight The next solution has put more material on the 
sup])C)it The direct connection between the seat and the bracket for placing the paddle is not present 
but it IS connected through the rear wheel bracket This arrangement reduces the stresses in this pait 
of the dc'sign The maximum weight solution is not having whole of the plate, instead it has placed the 
maternal more on the joint of the rear wheel braicket and the seat and the handle joint with the seat, to 
i('du( c the stresses 
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Weight 

Figuie 5 27 Hybiid proceduie to find five trade-ofF solutions for the bicycle frame design problem when 
weight IS not considered 





Figure 5 28 Five different designs of the bicycle frames when the weight of the frame is not accounted 
lor the loading 
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Weight 


Figme 5 29 Hybrid procedure to find nine trade-off solutions for the bicycle frame desigii problem when 
weight IS considered 




Figure 5 30 Five different designs of the bicycle frames when weight of the bicycle frame is also 
( onsidered 
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5.9 Closure 


Tlu' f hajiter piesents the results for the different multi-objective optimization test problems foi the cases 
u Ik' 11 the weight of the designed plate is considered while calculating the weight and when weight is not 
ajiplied to find the loadings The different problems are solved like — design of cross-section design oi 
( aiil ilevei plate design of simple supported plate with a number of loading cases, design of the hoistei 
platt' The design of the bicycle fiame is carried for a very high static load The hybiid method is used 
(oi solving diffeient pioblems The hybrid method uses combination of NSGA-II and hill climbing local 
s('ai( li 

The lesults obtained b^ the study are very interesting In all the design pioblems two confiictmg 
wc'ights aie used One is the minimization of the weight and other is the minimization of the deflection 
\aiying results fiom solutions ranging from the minimum weight (maximum deflection) solution to 
maximum weight (minimum deflection) solution, are obtained The results obtained give the designer 
innovative ideas foi the design of new shapes which have not been conceived so far A number of 
soliitioiis aic made available to the designer m a single simulation run The results for the real life 
liioblems like the design of the bicycle confirms the efficacy of the approach This makes the approach 
more' attractive and feasible to solve the real world design problems of varying complexity 
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Chapter 6 

Conclusion and Scope for future work 


riic pK'sont woik pioposes a new scheme of solving the shape optimization problems The design ol 
I 1k' shapes is earned out by considering a rectangular plate as the basic shape for the design Heie the 
sliajies aie lepiesented by the presence or absence of small connected elements, which aie repiesented 
bv the binaiy numbeis Thus the shape of any arbitrary structure can be represented by means of a 
sfiiiig ol biiiai\ numbers The shape obtained may have some disconnected regions To eliminate the 
disc oiiiKH ted legion the biggest cluster of the connected elements is found The shape lepresented bv 
t his iiK'thod may have some sharp corners To avoid that a smoothing method is devised and applied 
lieie to give the shape corresponding to the binary string 

The biiiaiy GA can handle the binary strings so this is chosen as the method of solving the shape 
dc'sigii pioblems Initially, a population of individuals is created randomly The binary individuals 
fiK' (ouvc'ited into the shapes according to the above described method The shapes which satisfy the 
gc'ouK'tiy ( oust 1 amts undergo finite element analysis, which give the maximum stress and maximum 
stiaiii values The solutions which satisfy the stress and strain constraints are considered feasible Thc' 
weight of these shapes is evaluated The solutions which are geometrically infeasible are assigned a 
voi y heavy penalty The shapes those are infeasible due to the stress or strain constraint violation are 
assigned an eiior value according to the violation of the constraint 

Foi single objective cases weight is taken as the fitness function and design is carried out for the 
minimization of the weight Now standard GA operators i e selection, crossover and mutation aie 
used to find new population of the individuals Here tournament selection with a constraint handling 
strategy is used for the selection procedure An innovative two dimensional crossover which respects 
the geometry of the individuals is devised and used This crossover is found to be efficient in these 
pioblems Single point mutation is used for mutation The solution is evolved through GA and get to 
a good solution from where, the local search is carried out on this best individual found by the GA 
When the problem is analysed for more than one conflicting objectives, the smgle objective GA is 
not suited as there are many optimal solutions The solution for these problems is a Pareto-optimal 
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s('t of solutions A specific multi-objective genetic algorithm — elitist non dominated genetic algorithm 
(NSGA-II) IS used to solve the multi-objective problems This algorithm has the capability of reaching 
1 1h' global Paioto-optimal fiont as well as obtaining the diversity on the front Heie two objectives aic 
iniuinu/at ion of tlio weight and the minimization of the deflection There are constiamts on the max- 
niinin clc'flection \alue and maximum stiess values The population undergoes the selection ciossovei 
.111(1 mutation Elitism is used to preseive the best non dominated solutions found so fai The Goldbeig s 
Milking scheme with constrained dominance is used to rank the population The selection is done on 
1 li(' basis of the lank and the diveisity on the Pareto-optimal front The two dimensional ciossovei and 
siiii])l(' mutation is used Each solution of the best solution set obtained by the NSGA-II undeigoes 
ilu' local seaich This results in a set of solutions, which may not be non-dominated Non-dommated 
solutions aie identified by a non-domination search on this set of solutions Fiom this final set of non 
dominated solutions, as many as requiied diverse solutions are selected using the clustciiug appioach 
The local seaich is done by flipping one bit at one time in an ordered manner Eveiy time the shapes 
<u(' c'xti acted fiom the new stung This is analysed and if the new solution is bettci than the pievious 
solutions, the change is accepted and if this is not, the change is rejected and the previous string is 
K'stou'd The local search technique uses a single objective strategy hence it poses no problem for the 
single objective pioblems For multi-objective problems a weighted sum strategy is used to convert the 
multiple objectives into single objective The weights are determined by using the fitness values of the 
solution Two stiategies are used for determining weights, one of the constant weights stiategy i e the 
wcnght vectoi is determined once, and second strategy is of continuous updating of the weights The 
optimum solution obtained after each local search is used to update the weight vectoi 

6.1 Conclusion 

A number of shape optimization problems axe solved using the hybrid approach The problems consid- 
('U'd aie design of cross-section, cantilever plate, simple supported plate, hoister design are considered 
A pioblem of design of bicycle frame is also used to show the application of the approach to the real 
woild problems The main results can be summarized as following - 

• The proposed hybrid approach is found to be very efficient to solve the shape optimization prob- 
lems 

• The solution depends on the representation scheme Proper representation is required to have the 
accuracy of the solution 

• The proposed elemental approach of representmg shapes is helpful in the finite element analysis, 
which also break the search domain into smaller elements Thus one computational step is reduced 
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• The use of elemental approach make it easier to find the solution for the cases when the properties 
of the material are varying over the search space 

• R)i the loal woild pioblems, where the optimal solution of the problem is not known it is safe to 
get lu'ai to optimal solution and then use of local search with a good solution to get the optimal 
solution This may piove an efficient approach to find the optimal solutions 

• It IS known that the classical local search methods has the potential to reach the optima if the 
initial guess is veiv good The same knowledge is used in this hybrid method to ensure the 
(oiivcigeiia' to the global optimal solution The genetic algorithm takes the solution to the basin 
of global optima and then local search starting from a good initial solution can take the solution 
to optima 

• Since the local seaich method reaches the optima very efficiently if it starts from a good solution 
so a consideiable amount of computational effort can be reduced by using the proposed approach 

• It IS found that since the solutions are only governed by the fitness values, many configurations 
fan be obtained foi the same fitness 

• The pi oposed method has the potential of finding the shapes which are obvious and results difficult 
to visuali/e can also be obtained without any difficulty 

• The proposed method can be very helpful in finding the new innovative solutions for the problem 
which may piove useful m conceptualizing new designs This feature of the considered approach 
make it oxtiemely useful for the designer 

• This appioach follows the evolution process to find the best possible shape thus mimicking the 
nature, hence the solution found are the best possible cases in most of the cases making the best 
use of the resources Hence this approach can be of interest m the wake of globalization, as the 
]ob of designei to make the best possible use of the resources is easily achieved 

• The appioach presented here give a more practical look towards the use of genetic algorithms 
particularly the multi-objective genetic algorithms m practice to solve the real world problems 

• For the multiple objective solutions it is possible to get the multiple solutions in one simulation 
lull 

• For the multi-objective cases the set of Pareto-optimal solution is obtained by the NSGA-II This 
is further improved by the local search and hence convergence to the global Pareto-optimal front 

IS ensured 
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1 1 he rippiOfUh proposed here shows a way of using the genetic algorithms for the real world 
jnohlenis It’s efficiency is well established through a number of test problems This call foi 
luiiiH'diatc' >is(> of the pioposed scheme to solve the real world optimization problems of stiuctuial 
t)])t iiui/.al 1011 

2 1 lu' <i])])ioa<li can be extended to the 3-diniensional cases by making use of thiee dimeiisioiial 

('h'liK'iits 

1 riu' pioposed appiodch uses a small number of elements to represent the shapes It will be 
ml ('lest mg to hnd the effect of increasing dimensionality of the problem 

1 rh(' piopos('d appioach uses first order polynomial for smoothing purposes It will be mtciestmg 
to liiul the applicability of higher order polynomials for this purpose 

a Tlu' optamal mesh generation techniques can be used to reduce the computational time 

() Oiii mt( lestmg task may be the parallelization of the process to make the method more piactual 
hv uHlucmg the computational time 

7 1 lu' solution of the bicycle design can be physically formulated and verified for the applicability 

ol t h(' (h'sign 

iS '1 1 h' (h'sign of the fragmented design and then coupling of all the designs to get the asseniblv 
(l('sigii mav be an interesting work 

<) Tlu' iiu't hod ( an be used to find the solutions when the properties of the material are not constant 
ov('i th(' (loniam. instead they are varying 

10 The' pioposed approach can be used for different optimization problems from other fields as well 

6.3 Closure 

l'h(' chajitc'i presents the brief sketch of the methodology used for representing the shapes and then 
appheataon of hybrid approach to solve the problems of shape optimization The results obtained from 
the' different test problems solved are used to draw certain conclusions The mam points axe presented in 
t,h(' conclusions The results show the applicability of the approach to the real world problems of shape 
dc'sign Tins approach can be used to solve the real world problems of the fields other than structural 
opinmzation The approach can also be extended to a more realistic case of 3-dimensional analysis 
The' futiiic^ scope of the work includes the applicability of the considered approach to the problems 
m otlu'i walks of life and also it asks for the physical manufacturing of the bicycle frame The extension 
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to i-diniciisioiidl cases and the use of fragmented design approach are also interesting fields for the 
flit iiK' woiks Use of parallel GA to find the solution of these problems will lead to the leduc tioii in the 
I oiiiputatioiial time This is lequiied to be pioved by the use of parallel stiuctiue of GA 
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